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Transverse momentum (pt) spectra of electrons from semileptonic weak decays of heavy flavor 
mesons in the range of 0.3 < < 9.0 GeV/c have been measured at midrapidity {\r]\ < 0.35) 
by the PHENIX experiment at the Relativistic Heavy Ion Collider in p+p and Au+Au colhsions 
at ^/»NN ~ 200 GeV. In addition, the azimuthal anisotropy parameter V2 has been measured for 
0.3 < py < 5.0 GeV/c in Au+Au collisions. The nuclear modification factor -Raa with respect to p+p 
collisions indicates substantial energy loss of heavy quarks in the produced medium. Comparisons 
of Raa and 112 are made to various model calculations. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

Numerous experimental results from the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Lab- 
oratory (BNL) have firmly established that the matter 
created in central Au+Au collisions at ^snn — 200 GeV 
cannot be explained by the common expectation of a 
weakly interacting gas of quarks and gluons at very high 
densities [H-Q- Strong suppression observed in measure- 
ments of the nuclear modification factor i?AuAu(PT) for 
light- fiavor hadrons shows that high-p^ scattered partons 
suffer a significant energy-loss in the medium [5|-[3]- The 
large magnitude of the differential elliptic flow parameter 
V2(pt) and its pr and mass dependencies measured in a 
limited phase space approximately agree with the theo- 
retical predictions based on an ideal hydrodynamic fluid 
[8l-[Tl| . In addition, the observation of universal scaling of 
V2 for various hadrons suggests that the flow pattern is es- 
tablished at the partonic level before hadronization [l^ . 
These experiments indicate that the medium is strongly 
interacting and exhibits hydrodynamic behavior. 

However, it has been pointed out that qualitative evi- 
dence for a near-perfect hydrodynamic fluid is not suffi- 
cient in terms of thermodynamic and transport concepts 
[T3I IT^ . The issue is how perfect is the near-perfect fiu- 
idity observed at RHIC. The validation of the perfect flu- 
idity requires the dimensionless ratio of shear viscosity 77 
to entropy density s to be small. 

Heavy quarks (charm and bottom) are important 
probes of the dense matter formed at RHIC. Because of 
their large masses, their dominant production mechanism 
is restricted to parton-parton collisions in the initial stage 
of the reaction. They can interact with the medium dif- 
ferently than light quarks and gluons due to their heavy 
mass. It was predicted that the energy loss of heavy 
quarks would be smaller than that of light quarks and 
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gluons due to suppression of small angle gluon radiation, 
called the "dead cone effect" [Ull^. 

Recently, heavy quark measurements along with mod- 
els for heavy quark interaction have opened new possi- 
bilities to investigate other interaction mechanisms such 
as coUisional energy-loss and in-medium fragmentation. 
In some models, the heavy quark diffusion coefficient 
D^'^ controls the extent to which the initial power-law 
Pt spectrum approaches the thermal spectrum and the 
extent to which the heavy quark will follow the under- 
lying flow of the medium. Simultaneous measurement of 
heavy-quark i?AuAu(PT) and V2{p?[) can provide an esti- 
mate of D^Q. 

Heavy quark measurements in p-\-p collisions serve as 
a testing ground of QCD. Because of their large mass, it 
is expected that next-to-leading order perturbative QCD 
(NLO pQCD) can describe the production cross section 
of charm and bottom at high energy, particularly at high 
Pt- At the Tevatron with ^/s = 1.9 TeV, bottom produc- 
tion is well described by NLO pQCD [l^. Charm pro- 
duction cross sections at high pT are found to be higher 
than the theory by « 50%, but are compatible within the 
theoretical uncertainties .18j. 

There are several ways to measure heavy quark pro- 
duction. The most direct method is to reconstruct D or 
B mesons from their decay products, such as D — > if tt or 
B — >■ J/ijjK. PHENIX is currently not capable of mea- 
suring the decay vertex of the heavy meson, which makes 
this method very difficult. The method which is em- 
ployed by this analysis is to measure single leptons from 
heavy flavor decay. Both the charm and bottom quarks 
have relatively large branching ratios 10%) to single 
electrons or single muons. Production yields and momen- 
tum distributions of the parent charm or bottom hadrons 
can be inferred from the invariant spectrum of single lep- 
tons from the decay. This method has an advantage in 
heavy-ion collisions at RHIC where the first method suf- 
fers from a very large combinatorial background due to 
the high multiplicity of the event. 

Single electrons in hadronic collisions were first ob- 
served in the early 1970's in -/s = 52.7 GeV p-\-p colli- 
sions at the CERN ISR J^] before the discovery of charm. 
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Subsequently, several experiments reported single elec- 
tron production at the ISR. In pp collisions at y^=630 
GeV, UAl [lO] measured bottom production via single 
muons and UA2 [HI reported the charm cross section 
from single electron measurements. At the Tevatron col- 
lider, the CDF and DO [1^ experiments measured 
bottom production via single muons or single electrons. 

At RHIC, PHENIX first measured charm production 
in Au-f-Au collisions at ,/s nn = 130 GeV via measure- 
ment of single electrons [2J| . Subsequently, PHENIX re- 
ported results of single electron measurements in p+p 
[25] and Au-|-Au collisions [26l - [28j at midrapidity from 
the 2002 data set. In Au-f Au collisions, the total yield 
of heavy-flavor decay electrons was found to scale with 
the number of binary nucleon-nucleon collisions [26| as 
expected for a point-like process. As already mentioned, 
the energy loss of heavy quarks was expected to be re- 
duced due to the dead-cone effect. Consequently, it 
was expected that high pT suppression and V2 of heavy- 
flavor electrons would be much weaker than those of light 
mesons. In contrast, a strong suppression of heavy-flavor 
electrons for 2 < < 5 GeV/c [28] and a non-zero 
electron V2 for px < 2 GeV/c were discovered. The 
STAR experiment has also measured the yield of elec- 
trons from heavy flavor decays at RHIC |29| . 

This article presents measurements of single electrons 
{Nf.+ + iVg- )/2 from semileptonic decays of heavy quarks 
(charm and bottom) at midrapidity (]77] < 0.35) in p + p 
and Au+Au collisions at \/snn=200 GeV. The data were 
collected in 2004 (Au-f Au), and 2005 {p + p). We extend 
the previous PHENIX analyses of electron measurements 
in p+p and Au-|-Au collisions at y^snn = 200 GeV to 
a broader p^ range: 0.3 < p^ < 9.0 GeV/c and with 
a much higher precision. Part of the results have been 
published in [33,|3l|. 

This article is organized as follows: Section HIl presents 
an overview of the PHENIX detector system related to 
the analysis. Section Hill presents the details of the data 
analysis. Section HVl shows the p^ distribution of the in- 
variant spectrum in the 2005 p+p and invariant yields 
and elliptic flow in 2004 Au-|-Au collisions. The results 
are compared with theoretical predictions and discussed 
in Sec.|Vl Finally, Sec. I VII gives a summary and conclu- 
sions of the analysis. 



II. PHENIX DETECTOR 

A detailed description of the complete PHENIX de- 
tector system can be found elsewhere [33 - l37| . Here we 
describe the parts of the detector system that are used 
in this analysis, namely, two global detectors and two 
central arm spectrometers. The global detectors are the 
beam-beam counter (BBC) and the zero-degree calorime- 
ter (ZDC). Each of the central arms covers a pseudora- 
pidity of \ri\ < 0.35 and an azimuthal angle of tt/2. They 
contain the drift chamber (DC) and the multiwire pro- 
portional pad chamber (PC) for charged particle track- 



ing, the ring-imaging Cerenkov detector (RICH) for elec- 
tron identification, and the electromagnetic calorimeter 
(EMCal) for energy measurement. Figure [1] shows the 
beam view of the PHENIX detector. 




Beam View East 



FIG. 1: (Color online) Beam view (at 2 = 0) of the PHENIX 
central arm detector in 2004 (Au-|-Au) and 2005 (p+p). The 
detectors used in the present analysis are the drift chamber 
(DC) and the multiwire proportional pad chambers (PCI, 
PC2, PCS) for charged particle tracking, the ring-imaging 
Cerenkov detector (RICH) for electron identification, and the 
electromagnetic calorimeters, which are Lead Glass (PbGl) 
and Lead Scintillator (PbSc), for energy measurement. 



A. Global detectors 

The BBC and the ZDC measure the collision time and 
the collision vertex position Zvtx along the beam. In 
Au-|-Au collisions, they are used to determine the cen- 
trality of the collision [s^]- They also provide the first 
level, minimum bias (MB) trigger information for beam- 
beam collisions. 

Two sets of BBCs are placed at ±1.44 m from the nom- 
inal interaction point along the beam line (one on the 
north side and one on the south). Each BBC consists of 
64 Cerenkov counter modules, arranged radially around 
the beam pipe. Each module is made of a 3 cm long 
quartz radiator and a photomultiplier tube (PMT). The 
PMTs can operate in a high magnetic field (^ 0.3 T) that 
lies mostly parallel to the beam axis. The BBC is used 
to estimate the number of charged particles in the pseu- 
dorapidity region 3.0 < \r]\ < 3.9. The vertex position 
resolution of the BBC is ~ 0.6 cm for Au-|-Au collisions 
and a few cm in p+p collisions. 

The ZDC measures the total energy carried by the for- 
ward neutrons produced along the beam direction. These 
neutrons are produced either by Coulomb dissociation of 
the beam particles or by evaporation from beam specta- 
tors. The ZDCs are placed at ±18 m from the interaction 
point along the beam line. The angular acceptance of 
each ZDC is \e\ < 2 mrad (]ry] > 6). Each ZDC consists 
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of three modules of two- interaction-length deep tungsten- 
quartz Cerenkov sampling calorimeters. The ener gy r es- 
olution of the ZDC is 5E/E - 2m/ (GeV) 



B. Central magnet 

The transverse momentum of each charged particle is 
determined by its bending curvature in the magnetic field 
provided by the PHENIX central magnet (CM) system 
[33| . The CM is energized by two pairs of concentric coils 
and provides an axial magnetic field parallel to the beam 
pipe. The coils can be run with the fields for the two coil 
sets adding (the "++" and " " configurations) or can- 
celing ("H — " configuration). During the Au-|-Au mea- 
surement in 2004 and the p + p measurement in 2005, 
the CM was operated in the "++" and " " config- 
urations. In the and " " configurations, the 

field component parallel to the beam axis has an approx- 
imately Gaussian dependence on the radial distance from 
the beam axis, dropping from 0.9 T at the center to 0.096 
T (0.048 T) at the inner (outer) radius of DCs. The total 
field integral is J B ■ dl ^ 1.15 T-m. 



C. Tracking detectors 

The DC and PC [s^l in the central arms measure 
charged particle trajectories in the azimuthal direction 
to determine the transverse momentum (pr) of each par- 
ticle. By combining the polar angle information from the 
innermost PC (PCI) with the px, the total momentum 
(p) is determined. 

The DC is positioned between 2.02 and 2.46 m in radial 
distance from the z axis for both the West and East arms. 
Each DC occupies a pseudorapidity of I77I < 0.35 and n/2 
in azimuth. Each DC volume consists of 20 sectors, each 
of which covers 4.5 degrees in azimuth. Each sector has 
six types of wire modules stacked radially, named XI, Ul, 
VI, X2, U2 and V2. Each module is further divided into 
4 drift cells in the (f> direction. A plane of sense wires is 
at the center of a drift cell, with 2 to 2.5 cm drift space. 
The X wires run parallel to the beam axis and measure 
the particle trajectory in the r-(f) plane. The U and V 
wires have a stereo angle of about 6.0 degrees relative to 
the X wires in order to measure the z-coordinate of the 
track. The single X wire resolution is ~ 150 fim. The 
intrinsic tracking efficiency of the X modules is greater 
than 99 %. Helium bags were placed between the beam 
pipe and the DCs to reduce the photon conversions and 
multiple scattering. 

The PC determine the spacepoints along the straight 
line particle trajectories outside the magnetic field. They 
are multiwire proportional chambers that form three sep- 
arate layers of the central tracking system. The first PC 
layer (PCI) is located between the DC and RICH at 2.47- 
2.52 m in radial distance from the interaction point, while 
the third layer (PCS) is located in front of EMCal occu- 



pying 4.91-4.98 m from the interaction point. The second 
layer (PC2) is placed behind RICH occupying 4. 15-4. 21m 
in radial distance in the West arm only. Position infor- 
mation from the PCI and the DC, along with the vertex 
position measured by the BBC, are used in the global 
track reconstruction to determine the polar angle of each 
charged track. 



D. Ring-imaging Cerenlcov detector (RICH) 

The RICH is a threshold-type gas Cerenkov counter 
and the primary detector used to identify electrons in 
PHENIX 35, 39]. It is located in the radial region of 
2.5-4. Im just outside PCI. The RICH in each central 
arm covers a pseudorapidity range of \ri\ < 0.35 and tt/2 
in azimuthal angle. Each volume contains mirror panels 
(0.53 % of a radiation length thick), forming two inter- 
secting spherical surfaces, with a total reflecting area of 
20 m^. The spherical mirrors focus Cerenkov light onto 
two arrays of 80((/>) x 16(z) = 1280 PMTs, each located 
on either side of the RICH entrance window as shown in 
Fig. [2] Each PMT has a magnetic shield that allows it 
to operate in a magnetic field up to 100 Gauss. CO2 gas 
at atmospheric pressure {n = 1.000410) was used as the 
Cerenkov radiator. The RICH has a Cerenkov threshold 
of 7 = 35, which corresponds to p > 20 MeV/c for an 
electron and p > 4.9 GeV/c for a pion. 



Most hadrons do not 
emit Cherenkov light 




Cherenl<ov ligtit 



FIG. 2: (Color online) Top view (at y — 0) of the ring-imaging 
Cerenkov detector (RICH) in the PHENIX East arm. 

The average number of hit PMTs per electron track 
is about 5, and the average number of photo-electrons 



detected is about 10. 
p + p. 



The c/tt separation is 



10" 



E. Electromagnetic calorimeter (EMCal) 

The EMCal [s^ consists of two arms, each covering a 
pseudorapidity region of \r]\ < 0.35 and tt/2 in azimuthal 
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angle. Each of the two arms consists of four rectangular 
sectors. Two adjacent sectors of the East arm are based 
on lead-glass (PbGl) calorimetry, whereas the remaining 
sectors are based on lead-scintillator (PbSc) calorimetry. 
The radial distance from the z axis is 5.10 m for PbSc 
and 5.50 m for PbGl. 

The PbSc is a Shashlik-type sampling calorimeter 
made of alternating tiles of lead and scintillator. It con- 
sists of 10.5 X 10.5 X 37 cm^ (18.2 Xq) rectangular mod- 
ules, each of which is constructed from alternating lay- 
ers of 1.5 mm thick lead and 4 mm polystyrene-based 
scintillator layers. The light from the scintillator is col- 
lected by optical fibers that run longitudinally through 
the module volume and are brought to four PMTs in 
the back end. This structure divides each module into 
four equal square cross section towers, from which the 
light is collected separately by the fibers. A PbSc sec- 
tor is subdivided into 18 (3 (0) x 6 (z)) super-modules 
made of 12 x 12 towers. Thus a PbSc sector consists of 
36((/)) x 72(z) = 2592 towers. The nominal energy reso- 
lution is 8EIE - 4.5 % ® (GeV) %. 

The PbGl consists of 4 x 4 x 40 cm^^ rectangular mod- 
ules, each assembled from 4 lead- glass crystals. Light 
emitted by particles depositing energy in the crystals 
is collected by a PMT at the back of each module. 
Twenty four of them (4 {(f) x 6 (z)) are combined 
together to form a PbGl-module. Six PbGl-modules 
(3 {(/>) X 2 (z)) comprise a PbGl-super-module and 32 
of them (4 (0) x 8 (z)) constitute a sector. This means 
a PbGl sector consists of 48(0) x 96(z) = 4608 PbGl 
counters. The PbGl h as a nomin al energy resolution of 
6E/E - 4.3 % ® 7.7//E7GeV) %. 

The EMCal served as the electron trigger device in the 
2005 p+p data taking. The trigger (PH) is fired whenever 
there is any hit which produces an energy of more than 
1.4 GeV in an EMCal trigger tile unit. The trigger tile 
consists of 4 X 4 non-overlapping towers. All trigger bits 
output from the tiles are summed in each super module 
online. Since there are 172 EMCal super-modules (108 
PbSc + 64 PbGl), the PH trigger has 172 output chan- 
nels. 95 % of them were live in the 2005 p-^p data. All 
the output bit information is summed and used for the 
event-trigger decision. The trigger efficiency is evaluated 
in section UlI Gl 



III. ANALYSIS 

A. Data set and event selection (Au+Au) 

The data for Au-|-Au collisions at ^snn = 200 GeV 
were taken during the 2004 run of RHIC. The minimum 
bias (MB) trigger was defined as: 



MB = (BBC > 2)&&(|z| < 38 cm) 



(1) 



trigger also requires at least one hit in one of the ZDCs. 
The vertex position along the beam line is determined 
online from the timing difference between the two BBCs. 
The MB trigger efficiency for inelastic Au-|-Au collisions 
is evaluated as 92.2j:3;o%. 

In the offline analysis, a tighter vertex cut (|bbcz| < 
20 cm) is required to eliminate conversion electrons from 
forward materials. Here bbcz is the vertex position de- 
termined by BBC in the offline analysis. Any abnormal 
data due to unusual beam or detector conditions are re- 
moved. Table |T] summarizes the number of MB events 
used for the analysis of electron spectra. For the V2 anal- 
ysis, only those events with a good reaction plane cali- 
bration were used. The number of MB events used for 
the V2 analysis was 7.1 x 10^ events. 

The MB events are divided into five centrality classes 
of 0-10 %, 10 -20%, 20-40%, 40-60%, and 60-92%. See 
section UlI CI for the centrality determination. 

There are two groups of runs: normal runs and con- 
verter runs. In the converter runs, an additional photon 
converter was inserted around the beam pipe. See sec- 
tion nil Jl for the converter runs and their use for the 
analysis. 



TABLE L Number of events for each centrality class for runs 
with and without additional converter material installed. 



Centrality class 



Normal runs 



Converter runs 



0-92% 
0-10% 
10-20% 
20-40% 
40-60% 
60-92% 



7.48 X 10** 

8.06 X lO'^ 

8.07 x lO'^ 
1.61 x 10* 
1.61 x 10* 
2.65 x 10* 



5.79 x 10^ 
6.31 X 10^ 
6.29 X 10^ 
1.25 X 10^ 
1.24 X 10^ 
2.04 X 10^ 



Here, BBC > 2 means that at least two hits are re- 
quired in both the North and South BBCs. The offline 



B. Data set and event selection (p + p) 

The data for p-\-p collisions at -^i = 200 GeV were col- 
lected during the polarized p + p collisions in 2005. The 
proton beams had approximately 50% longitudinal po- 
larization with alternating spin orientations in successive 
bunches. The polarization of the protons has negligible 
effect on the cross section. 

Two data sets are used for the electron analysis: (1) 
the minimum bias (MB) data set and (2) a "photon" 
(PH) trigger data set. The MB trigger for p + p required 
at least one hit in both the North and South BBC de- 
tectors in coincidence with the beam bunch crossing and 
the event vertex position within |bbcz| < 30 cm from 
the nominal collision point along the beam axis. The PH 
trigger required a minimum energy deposit of 1.4 GeV in 
an overlapping tile of 4 x 4 EMCal towers in coincidence 
with the MB trigger. The PH trigger had nearly 100% 
trigger efficiency for electrons with px above 2 GeV/c in 
the active trigger tiles. 
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The MB trigger cross section is (Tbbc = 23.0 ±2.2 mb. 
Since only ~ 50% of inelastic p + p collisions satisfy the 
MB trigger condition, only a fraction of the inclusive elec- 
tron production events were triggered by the PH trigger. 
This fraction is assumed to be momentum and process 
independent, and is determined to be ebias = 0.79 ± 0.02 
[4^ from the yield ratio of high pT tt'^'s with and without 
the BBC trigger. The value of ebias is slightly higher than 
that of 2002 (0.75) [H, due to a lower BBC threshold. 

After the selection of good runs and the vertex cut, an 
integrated luminosity (£) of 45 nb~^ in the MB data set 
and 1.57 pb~^ in the PH data set are used for the analy- 
sis. During a part of the 2005 p + p run, the same photon 
converter that was used in the 2004 Au+Au run was in- 
serted around the beam pipe. The integrated luminosity 
in the converter run period is approximately 7% of the 
total luminosity. 
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FIG. 3: (Color online) The correlation between fractional 
BBC charge Q/Qmax and fractional of ZDC energy E/Emax. 
The sections from right to left correspond to centrality classes 
0-10 %, 10-20 %, 20-40 %, 40-60 %, and 60-92 %. 



C. Centrality (Au+Au) 

The centrality of each Au+Au collision is determined 
by the measurements of BBC charge and ZDC energy. 
Figure [3] shows the correlation between BBC total charge 
and ZDC total energy for real data, each of which is nor- 
malized by the maximum value. The amount of BBC 
charge is proportional to the particle multiplicity and is 
correlated to the overlapping area of two colliding nuclei. 
In contrast, the total ZDC energy is smaller in peripheral 
collisions than it is in central collisions. In the most pe- 
ripheral collisions, a large fraction of the nuclei become 
spectators. The number of free neutrons reaching the 
ZDC becomes small, due to the formation of deuterons 
or heavier charged nuclei from the spectator neutrons 
and protons, which are swept by the magnets. The cen- 
trality determination from the BBC and ZDC energy is 
schematically shown in Fig. [31 Each centrality class con- 
tains approximately the same number of events . The 
most peripheral class is normalized to be 92%. 

The relation between centrality, the number of binary 
collisions (A'coii), and the number of participants ( N„„rt) 
is obtained from a Glauber Monte Carlo simulation [4l| | . 
The systematic uncertainties in A'coii and A^part are cal- 
culated from the uncertainty in the Glauber parame- 
ters, the centrality determination, and BBC and ZDC 
responses. Table HIl summarizes the average A'coii, A'pait, 
the associated nuclear overlap function (TauAu) and their 
systematic errors for each centrality class. 



D. Track Reconstruction 

Charged particles are reconstructed by the DC and 
PCs. Figure S] shows a trajectory up to PCS and asso- 
ciated kinematic parameters. Assuming that the track 
originates on the beam axis, the intersections of the tra- 
jectory with each detector plane are uniquely determined 
by four initial kinematic parameters: Zvtx, 9, (f'v, and the 



TABLE II: iVpart, A''coii and TauAu [iH by Glauber calculation 
for Au+Au collisions at a/snn = 200 GeV. 



Centrality class (A^part) (syst) {N^m) (syst) (rAuAu)(MB-') (syst) 



0-92% 


109.1 (4.1) 


257.8 (25.4) 


6.14 (0.45) 


0-10% 


325.2 (3.3) 


955.4 (93.6) 


22.75 (1.56) 


10-20% 


234.6 (4.7) 


602.6 (59.3) 


14.35 (1.00) 


20-40% 


140.4 (4.9) 


296.8 (31.1) 


7.07 (0.58) 


40-60% 


59.95 (3.6) 


90.70 (11.8) 


2.15 (0.26) 


60-92% 


14.50 (2.5) 


14.50 (4.00) 


0.35 (0.10) 



(a) Beam View 




(b) Side View 




FIG. 4: (Color online) A charged particle trajectory and the 
kinematic parameters are shown in (a) beam view, (b) side 
view, and (c) top view of PHENIX central region. 
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transverse momentum px- Here, z^tx is the interaction 
vertex along the z-axis, 6 is the angle between the ini- 
tial direction of the particle track and z-axis, and 0v is 
the initial azimuthal angle of the trajectory. These initial 
parameters are reconstructed frommeasured variables: a, 
P, Zpad, and (/)DC- ot is the angle between the projection 
of the trajectory in the x — y plane and the radial direc- 
tion, at the intersection point of the trajectory with the 
circle of DC reference radius (i?DC = 2.2 m). Zpad is the 
intersection point of the trajectory with PCI surface ra- 
dius (i?pci — 2.48 m). (3 is obtained by considering the 
plane which includes the z-axis and Zpad- It is defined 
as the angle between the projection of the trajectory to 
the plane and the line connecting the vertex to the posi- 
tion of the pad chamber hit. 0dc is the azimuthal angle 
of intersection point of the trajectory with the circle of 
radius with i?DC- 

The tracking process starts by collecting hits in DC 
XI and X2 wires. The hits are projected to the x — y 
plane at z = 0. Then, a Hough transform is applied to 
find the sets of hit points based on the assumption of 
the track having a straight line trajectory inside the DC 
volume [4^ |43| . The technique is performed using all 
possible XI and X2 hit combinations taking a and (/)dc 
as the parameters in Hough space. After this process, 
the direction of the found track line is specified by a and 
^DC- The hits in the U and V wires and PCI hits are 
associated with the track in order to obtain Zpad and /3. 
After the pattern recognition and track reconstruction, 
the initial momentum vector of the track at the collision 
vertex is calculated. A look-up-table technique is used 
for fast processing [H, . After the initial kinematic 
parameters are obtained, each reconstructed track is as- 
sociated with hit information from the outer detectors 
(PC2, PCS, RICH, and EMCal). 

E. Electron Identification 

Electron candidates are selected from reconstructed 
charged tracks based on information from the RICH and 
EMCal. For each of the reconstructed tracks, the follow- 
ing variables are calculated and used for electron identi- 
fication (elD); 

• quality - Tracks are required to have hits in both 
the XI and X2 sections of the DC and be uniquely 
associated with hits in the U or V stereo wires. At 
least one matching PCI hit is also required. 

• emcsdphi_e - Displacement in </) of the electron hit 
position of the associated EMCal cluster from the 
projected position of the track in units of standard 
deviations. For example, emcsdphi_e < 2 means 
that the position of the associated EMCal cluster in 
(j) is within 2 cr of the projected track position. This 
variable was calibrated specifically for electrons. 

• emcsdz_e - Same as emcsdphi_e, but for the z 

coordinate. 



TABLE III: Electron ID cuts used in Au+Au analysis 

elD cuts for Au+Au 
\/emcsdphi_e^ + emcsdz_e^ < 2.0 
nO > 2 
nl > 3 if pT < 5.0 GeV/c 
nl > 5 if pT > 5.0 GeV/c 
disp < 5.0 
chi2/npe0 < 10.0 
prob > 0.01 if pT < 5.0 GeV/c 
prob > 0.2 if PT > 5.0 GeV/c 
dep > -2.0 



• prob - The probability that the associated EMCal 
cluster is an electromagnetic shower. This variable 
is calculated from the value between the actual 
tower energy distribution of the cluster and the ex- 
pected distribution for an electromagnetic shower. 
For example, a cut, prob > 0.01, has 99% efficiency 
for a photon or electron shower, while it rejects a 
large fraction of hadrons. 

• nO - Number of hit RICH PMTs in an annulus 
region with inner radius of 3.4 cm and outer ra- 
dius of 8.4 cm around the the track projection on 
the RICH. The expected radius of a Cerenkov ring 
emitted by an electron is 5.9 cm. 

• nl - The number of hit RICH PMTs within 11 cm 
around the projection point of the track. 

• chi2/npe0 - A x^-like shape variable of the RICH 
ring associated with the track. 

• disp - A variable representing the displacement of 
the RICH ring center from the projected track po- 
sition. Units are cm. 

• dep - A variable of energy momentum match- 
ing. This variable is calculated as dep = {E/p — 
l)/aE/p, where E is the energy measured by EM- 
Cal; p is the momentum of the track; and <Je/p is 
the standard deviation of Gaussian- like E/p distri- 
bution. E/p is less than 1 for hadrons since hadrons 
do not deposit their full energy in the calorimeter. 
a E/p depends on the momentum of the electron. 

Tables IIIII and IIVI summarize the elD cuts used for 
the Au+Au and p + p analyses, respectively. For tracks 
with pt below 4.8 GeV/c, the RICH is fired only by elec- 
trons. Thus the RICH is the primary means of electron 
detection in this pr region. The E/p cuts remove re- 
maining background. The main cause of electron mis- 
identification below 5 GeV/c is that a charged hadron 
track is accidentally associated with hits in the RICH 
detector. Since the particle multiplicity in p + p is very 
small, typically few tracks in each of the central arm spec- 
trometers, the probability of such an accidental overlap 
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TABLE IV: Electron ID cuts used in p + p analysis 

elD cuts for p + p 
emcsdphi_ej < 4 
emcsdz_ej < 4 
nO > 2 
nl > 2 
prob > 0.01 
0.50 < E/p < 1.3(0.2 < p < 0.3 GeV/c) 
0.55 < E/p < 1.3(0.3 < p < 0.4 GeV/c) 
0.60 < E/p < 1.3(0.4 < p < 0.6 GeV/c) 
0.65 < E/p < 1.3(0.6 < p < 0.8 GeV/c) 
0.70 < E/p < 1.3(0.8 < p < 5.0 GeV/c) 



is much smaller. This permits us to use looser electron 
cuts in the p + p analysis. Note that although hadrons 
could fire the RICH above 4.8 GeV/c, the number of 
photoelectrons is small for hadrons below 5 GeV/c, and 
such hadrons are easily distinguishable from electrons. 
Thus we tighten our cuts above 5 GeV/c without worry- 
ing about extra hadron contamination. 

In p + p collisions the hadron contamination after the 
elD cuts is very small. The contamination level is esti- 
mated by reversing the prob cuts to enhance the hadron 
background. The estimated hadron contamination is 3% 
at pt = 0.3 GeV/c and less than 1% for 0.8 < pr < 5 
GeV/c with elD efficiency of approximately 90%. In 
Au-|-Au, there still remains hadron background due to 
accidental overlap between a track and RICH hits. This 
background is estimated and subtracted as described be- 
low. 

Figure [5] (a) and (b) show E/p distributions for two 
typical PT ranges in Au-|-Au collisions. All cuts in ta- 
ble |TTT] are applied except for the dep cut since it is just a 
PT-dependent cut in E/p. The square points are the esti- 
mated distributions of the remaining hadron background 
which are randomly associated with a ring in the RICH. 
This estimation is performed by swapping the north and 
south sides of the RICH in the same event in software, 
and reconstructing the track matching to RICH. That is, 
DC tracks from the south are matched with RICH hits 
in the north, and vice-versa. Since the north and the 
south sides of the RICH are identical and there was only 
~ 1 percent of dead channels in the RICH PMTs, this 
method gives a proper statistical estimate of the random 
liadron background in the electron sample. 

After subtraction of the random hadron associations 
(shown as squares in Fig. [5] (a) and (b)), an additional 
low E/p tail remains in the distribution at high px, as 
shown in Fig. [S] (b) . This tail is due to electrons from 
kaon decay (ifes) and photon conversions that occur far 
from the collision vertex. These background electrons are 
reconstructed with a momentum higher than the actual 
momentum of the electrons, and as such have a low E/p. 
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FIG. 5: (Color online) E/p distributions for various pT ranges, 
(a) and (b) are for Au-|-Au collisions, (c) is for p -I- p colli- 
sions, where the black stars are the data points and the solid 
circles show the GEANT [4^ simulation. The triangles and 
squares show the contributions from simulation of tt" and Kes, 
respectively. See text for more details. 
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A full GEANT simulation of the PHENIX detector 
was performed to determine the background from i^es de- 
cays and from photon conversions. Figure EJc) shows the 
E/p distribution in p + p for 0.8 < pr < 1.0 GeV/c com- 
pared with the GEANT simulation. The black stars are 
the data points and the circles show the simulation. The 
triangles and squares show the contributions from simu- 
lation of tt" and K^s, respectively. The simulated tt" and 
Ke3 decays went through the same offline analysis chain 
as in the real data, and identical electron identification 
cuts were applied. In the p+p data, the contamination of 
hadron mis-identification is negligible. In the tt'^ simula- 
tion, electrons are mainly produced either by the Dalitz 
decay or by photon conversion in the beam pipe. These 
conversion electrons produce the peak around E/p = 1. 
Note that in Fig. ^c) the data is above the simulation 
since the data also contain nonphotonic electrons. A 
small number of conversions in the helium bag contribute 
to the tail in the low E/p region. The main cause of the 
low E/p tail is Kg^ decay, shown as squares. Since the 
simulation reproduces the low E/p tail of the data very 
well, we conclude that the Ke3 background under the 
Gaussian peak, which will remain after the dep or E/p 
cuts, can be determined from the simulation. 

There is no adjusted normalization in this simula- 
tion/data comparison. In the simulation, we assume 
da/dy{pp tt") = 46.0 mb, da/dy{pp K^) = 
da/dy{pp — !■ K^) — 4.0 mb at midrapidity. The momen- 
tum distribution of tt" and kaons are based on PHENIX 
measurements in p + p. 

For pt above 4.8 GeV/c, charged pions begin to radi- 
ate Cherenkov light in the RICH. Thus tighter electron 
selection cuts are applied to reject pion background. In 
both the Au+Au and p + p analyses, tighter cuts of nl 
> 5 and prob > 0.2 are added. In the p + p analysis, 
the E/p cut is tightened to 0.8 < E/p < 1.3. With these 
cuts, the electron measurement is extended to 9 GeV/c 
in Pt- The tighter cuts are applied to tracks above 5 
GeV/c, as the Gherenkov radiation is very weak between 
4.8 and 5 Gev/c, and the data were binned in multiples 
of 0.5 GeV/c. 

The remaining hadron background with the tighter 
cuts is studied using the shape of the E/p distribution. 
Here we rely on the fact that the distribution of the 
prob and E/p variables are roughly independent of pT 
for hadrons at high pr, and that a cut on prob<0.01 
eliminates the vast majority of electrons. First, we ob- 
tain a sample of hadrons in the pr range of 1-4 GeV/c 
by imposing a veto on the RIGH. The hadron sample 
is then divided into two samples, one with prob>0.01 
and the other with prob<0.01. The ratio of these two 
hadron samples is taken. The upper plots in Fig. [S] show 
the two hadron samples and their ratio. The E/p distri- 
butions for Pt above 5 GeV/c are then estimated from 
the data with cuts identical to those described above, ex- 
cept with the prob cut reversed to prob<0.01. These 
E/p distributions are then divided by the ratio shown in 
the upper right plot of Fig. [B] to obtain an estimate of 



the E/p distributions of hadrons passing the prob>0.01 
cut. 

The lower plots of Fig. [5] show the estimate for the 
hadronic background along with the total E/p distri- 
bution. This procedure actually gives an overestima- 
tion of the background, as some electrons do pass the 
prob < 0.01 cut. This can be seen in the lower left plot 
as a peak in the hadronic sample. From these plots, we 
determine that any hadronic background below 8 GeV/c 
is negligible. Between 8-9 GeV/c in pT, we estimate a 
background of 20%, with an uncertainty of 10%. The 
elD efficiency of the tighter cuts is pt independent. In 
p+p, it is determined to be 57% of that for pT < 5 GeV/c 
by applying the same tighter cuts for pT < 5 GeV/c. 



F. Acceptance correction 

The acceptance correction in this analysis covers the 
following three components: the geometrical acceptance 
correction (e^®°), the elD efficiency (e°^^) and the recon- 
struction efficiency (e''°™). e^'^" accounts for the fraction 
of electrons that do not hit any detectors due to the finite 
solid angle covered by the detectors. e°^^ is the correc- 
tion factor for signal loss by elD cuts. e''°'^° takes into 
account that the measured electron spectrum in the de- 
tector is different from the real spectrum due to detector 
responses and the track reconstruction method. All the 
values are computed at once by the full detector sim- 
ulation for single events. Single and e~ events 
are generated uniformly in phase space {0 < p^ < 15 
GeV/c, \y\ < 0.5, and < < 2tt). They are processed 
by the fuh GEANT simulation program of the PHENIX 
detector [46j. The output simulation data files are pro- 
cessed by the event reconstruction chain of PHENIX. 
The acceptance correction factor (e^p is calculated as 
follows: 



^ICCO 



dN°''^/dp^ 
dNl^'/dp^ 



(2) 



The p^ distributions of input and output electron 
yields are dN^^/dp^ and dNf^^/dp^, respectively. The 
same elD cuts and fiducial cuts used in real data pro- 
cessing are applied to the output. Since the input elec- 
trons are generated flat in we must take into account 
that the smearing of the spectrum due to imperfect mo- 
mentum reconstruction affects a steeply falling spectra 
differently than our simulated spectra. Each track re- 
constructed in the simulation is weighted by the p^ of 
the corresponding input electron. 

Figure [7] shows three curves of as a function of 
Pt for (iVe+ -I- Ng-)/2. The top, middle and bottom 
curves are calculated with each set of elD cuts applied 
for 0.3 <PT < 5.0, 5.0 < pt < 7.0, and 7.0 < pt < 9.0 
GeV/c respectively. 
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Hadrons from RICH veto with 1 GeV/c < < 4 GeV/c 
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FIG. 6: (Color online) E/p distributions for (upper) hadron samples from a RICH veto cut which are used in the contamination 
study and (lower) the indicated high-pr ranges, shown along with the estimation of the hadronic contamination. 



G. Trigger efficiency (p + p) 



analysis. 



The trigger efficiency of the PH trigger in p+p colli- 
sions is determined from the MB data set as the fraction 
of electrons that satisfies the PH trigger. Since the trigger 
condition determined by the online level- 1 trigger proces- 
sor was recorded, we can apply the same PH trigger to 
the events recorded in the MB trigger in the offline anal- 
ysis. We require that the PH trigger bits of the event are 
set, and that the trigger tile that fires the PH trigger is 
hit by the electron candidate track selected by the offline 



Figure |S] shows the trigger efficiency of the PH trigger 
thus determined. The effective trigger threshold is about 
1.4 GeV, and the efficiency saturates above 2 GeV. The 
trigger efficiency at the plateau is « 86%, consistent with 
the fraction of active trigger tiles. The curve is a Fermi- 
like function that is fitted to the data. The fitted function 
is used as the trigger efficiency in the analysis. 
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FIG. 7: (Color online) e^l for (iV^+ + A^^- )/2 as a function of 
Pt- The top, middle and bottom curves are calculated with 
elD cuts applied for 0.3 < pr < 5.0, 5.0 < pr < 7.0, and 
7.0 < Pt < 9.0 GeV/c respectively. 
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FIG. 8: (Color online) Efficiency of the PH trigger in p + p 
determined from the minimum-bias data set. 



H. Occupancy correction (Au+Au) 

In addition to the efficiency for a single electron passing 
through the detector to pass the fiducial and elD cuts due 
to detector dead area and the efficiencies described above, 
there is also a finite efficiency loss for particle detection 
due to the presence of other particles nearby. To get a 
quantitative understanding of the multiplicity-dependent 
efficiency loss, we embed simulated single electrons and 
positrons into data files containing detector hits from real 
collisions. The same simulated particles are used in this 
method as those used to estimate the single particle ef- 
ficiency loss. The simulated particles are embedded into 
events such that the z vertex of the simulated particle 
and the event are similar. The simulated are run 
through the GEANT simulator of PHENIX, and the hits 
are added to the data files containing hits from a real 
Au-f Au event. Next, these new files containing the em- 
bedded are run through the entire reconstruction soft- 



ware to produce track candidates containing the variables 
upon which we make identification cuts. We then define 
the embedding efficiency as 



^embcd 



^ reconstructed from embedded data 
# reconstructed from single track data 



(3) 

where a reconstructed particle from embedded data has 
most of its detector hits associated with hits from the 
simulated particle. 

As a systematic cross-check on this method for deter- 
mining the multiplicity-dependent efficiency loss, a data- 
driven method was employed. The general strategy of 
this data-driven method is to select a very pure sample 
of electron conversion pairs from the data with an in- 
variant mass cut. We apply tight elD cuts to one of the 
tracks in a given pair to increase the chance that the 
pair really is an electron pair. Then, we measure the ef- 
ficiency loss as a function of collision centrality of tighter 
cuts (cuts actually used in the analysis) relative to that 
of the loose cuts. We then still need simulation to deter- 
mine the multiplicity-dependent loss due to the loose elD 
cuts, but simulation is more reliable when we use loose 
cuts. In order to get a pure sample of electrons, we iso- 
late pairs from photon-conversions in the beampipe 
and from tt'^ Dalitz decays. We assume, reasonably, that 
the multiplicity-dependent efficiency loss for these elec- 
trons is the same as for all electrons. 

30000 
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FIG. 9: Invariant mass distribution of e pairs from min-bias 
Au+Au data used for the multiplicity-dependent efficiency 
loss estimate. 

Figure [9] shows the invariant mass distribution of 
pairs from the data. The peak at ^ 20 MeV/c^ corre- 
sponds to the apparent mass of electrons from photon 
conversions in the beampipe. The peak at ^ 5 MeV/c^ 
corresponds to from tt" Dalitz decays. As mentioned 
above, we would like to have electrons identified by using 
both loose cuts and tight cuts. However, if loose cuts 
are applied to both the electron and the positron in the 
pair, then the sample will contain contamination from 
hadrons. In order to make sure that the pair is really 
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an e pair, we apply tight elD cuts to one of the parti- 
cles in the pair. We restrict ourselves to pairs with one 
of the particles in the acceptance of the Time of Flight 
detector ^3§\ so that we can use the TOF for electron 
identification. 

The cuts used for the analysis for pr < 5 GeV/c are 
then applied to the other particle in the pair, and the 
relative efficiency loss is calculated. As described in the 
previous section, for px > 5 GeV/c we obtain the effi- 
ciency by the relative electron yield between the samples 
with the cuts used for the two Pt ranges. This rela- 
tive yield includes multiplicity-dependent effects, so we 
only explicitly calculate the multiplicity-dependent effi- 
ciency loss for the cuts used for pT < 5 GeV/c. Some 
of the relative efficiency loss using the method just de- 
scribed is due to the single particle efficiency. To obtain 
an estimate of the multiplicity-dependent component of 
this relative loss, we divide the relative loss by the loss 
for peripheral collisions (60%-93% centrality), since the 
multiplicity-dependent loss in peripheral collisions is very 
low. So we have defined the multiplicity-dependent effi- 
ciency loss from the data-driven method for centrality C 
as 

data.C 

data,C _ ^cID fA\ 

^mult ~ data,60 93% ^ ' 

^cID 

with 

data.C ^tight cuts 

^elD — ~C ■ ^"^f 

loose cuts 

Ideally we would measure the above efficiency with no 
bias caused by making the initial loose cuts. We correct 
the multiplicity-dependent efficiency loss by using the 
efficiency from the embedding simulation for the loose 
cuts, e^j^gg'^''^ to obtain the more accurate estimate of 
the multiplicity-dependent loss from the otherwise data- 
driven method. 

TABLE V: Multiplicity-dependent efficiency loss calculated 
by two methods. £<=™t'<=d,c ^j^g efficiency loss calculated from 
the embedding simulation. e[^^^°'*''^ is the same as above but 
with the loose cuts. 



Centrality 


^cmbcd,C 


data,C embed, C 
^mult ^loosc 


00-93% 


0.852 


0.882 


0-10% 


0.771 


0.769 


10-20% 


0.835 


0.856 


20-40% 


0.900 


0.924 


40-60% 


0.952 


0.977 


60-93% 


0.982 


0.997 



Table |V] displays the various embedding efficiencies de- 
scribed above. The difference between the left and right 
columns gives an idea of the systematic error involved in 
this estimation. The use of the sample of conversion and 
Dalitz pairs will be discussed further in the section 
pertaining to systematic error analysis. 



I. Electron background cocktail 

The inclusive electron spectra consists primarily of 
three components: (1) "nonphotonic" electrons from 
heavy-flavor decays, (2) "photonic" background from 
Dalitz decays of light neutral mesons and photon conver- 
sions (mainly in the beam pipe), and (3) "nonphotonic" 
background from K — env (isTes), and dielectron decays 
of vector mesons. 

At high electron px, Drell-Yan processes also account 
for a small but non-negligible contribution to the elec- 
tron spectrum. Also, heavy quarkonia decays contribute 
at high electron px, which, although arguably can be in- 
cluded in electrons from "heavy-flavor" decays, should 
be distinguished from open heavy flavor decays for the 
purpose of interpretating the measurement. 

The photonic background is much larger than the 
nonphotonic background. The signal of electrons from 
heavy-flavor decays is small compared to the photonic 
background at low px {S/B < 0.2 for px < 0.5 GeV/c) 
but rises with increasing px {S/B > 1 for > 2 GcV/c). 
In order to extract the heavy-flavor signal, the back- 
ground has to be subtracted from the inclusive electron 
spectra. One technique to accomplish this task is the so- 
called "cocktail subtraction" method described in detail 
here. 

A cocktail of electron spectra from background sources 
is calculated using a Monte Carlo event generator of 
hadron decays and then subtracted from the inclusive 
electron spectra. This technique requires that the phase 
space distributions of the relevant background sources be 
well known. The PHENIX measurements of the relevant 
electron sources are precise enough to constrain the back- 
ground within a systematic uncertainty better than 15 % 
for all px- This uncertainty is of the same order as the 
signal to background ratio at the lowest px and, there- 
fore, it is not sufficiently small to extract the heavy-flavor 
signal via the cocktail subtraction over the full px range. 
Hence, at low px, a complementary technique to sub- 
tract the background, the so-called "converter subtrac- 
tion" method, is used to extend the heavy-flavor mea- 
surement to the lowest px with good precision. Conse- 
quently, the converter subtraction is the key to extracting 
the total heavy-flavor yield or cross section since most 
of the electrons from heavy-flavor decays have low px- 
However, towards high px, e.g. for px > 2 GeV/c where 
the converter subtraction starts to suffer from a lack of 
statistical precision, it is beneflcial to apply the cocktail 
subtraction since the signal to background ratio is large, 
statistics is irrelevant due to the Monte Carlo nature of 
the cocktail subtraction, and the cocktail input is known 
with small systematic uncertainties as discussed in the 
following. 

The most important background source, except for the 
highest electron px where contributions from direct radi- 
ation become dominant (see below) , comes from the neu- 
tral pion. The contribution of tt" decays to the photonic 
background is twofold. First, the Dalitz decay of neutral 
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FIG. 10: (Color online) (upper) Invariant differential cross 
section of charged pions (squares at low pr) and neutral pions 
(circles) for p + p collisions, together with a fit according to 
Eq. [S] (lower) Ratio of the data to the fit. 



pions (tt" — )■ e+e~7) is a primary source of electrons froni 
the collision vertex. Second, the conversion of photons 
from the decay 7r° — 77 in material in the PHENIX cen- 
tral arm aperture (mainly in the beam pipe) gives rise 
to a secondary source of electrons originating away from 
the original collision vertex. It is crucial to note that the 
contribution from photon conversions is smaller than the 
contribution from Dalitz decays. This is due to the care- 
fully minimized material budget in the PHENIX central 
arms. Apart from the beam pipe, which is made out of 
Beryllium and contributes less than 0.3 % of a radiation 
length to the material budget, helium bags constitute the 
only material between the beam pipe and the tracking 
and electron identification detectors in PHENIX. As was 
verified in a full GEANT simulation of vr" decays, the 
ratio of electrons from the conversion of photons from 
7^0 'j'y decays to electrons from 7r° Dalitz decays is 

0. 403 with a systematic uncertainty of about 10 %. This 
ratio is independent of pt in the range relevant here, 

1. e. for Pt > 1 GeV/c. For heavier mesons this ratio is 
rescaled in the cocktail to properly account for the fact 
that the branching ratio for the Dalitz decay relative to 
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FIG. 11: (Color online) (upper) Invariant multiplicity of 
charged pions (squares at low pr) and neutral pions for 
Au+Au collisions from 2002 (triangles) and 2004 (circles) to- 
gether with a fit according to Eq. [S] (lower) Ratio of the data 
to the fit to the 2004 data. 



the 77 decay grows slightly with increasing parent meson 
mass. Consequently, once the phase space distributions 
of tt'^ and heavier mesons are available it is straight for- 
ward to determine both the Dalitz decay contribution 
to the background electron spectrum as well as the cor- 
responding contribution from the conversion of photons 
from the same parent mesons. 

The phase space distributions of tt'^ are obtained via 
simultaneous fits to the charged and neutral pion spec- 
tra. This approach is only valid under the assumption 
that the invariant 7r° spectra and the averaged charged 
pion spectra (7r+ + Tr~)/2 are the same. While this as- 
sumption in general is well justified, at low px, i.e. for 
Pt < 1 GeV/c, the decay of ij mesons into three 7r° cre- 
ates a tiny charge asymmetry. However, according to a 
PYTHIA [13] calculation and consistent with data, this 
asymmetry is only about 2% and can be safely ignored 
in this context. 

Figure [TU] shows the comparison of the neutral and 
charge-averaged invariant differential cross sections of pi- 
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ons in p+p collisions at y/s = 200 GeV in comparison 
with a simultaneous fit to the data with a modified Hage- 
dorn parametrization: 



d^p {exp{-apT - bp'^) + pr/poY 

where a, 6, c, and n are fit parameters. Both an 
absolute comparison as well as the ratio of the data to 
the fit are shown to demonstrate the excellent quality of 
the parametrization. 

For Au+Au collisions the 7r° invariant differential mul- 
tiplicity distributions are obtained by equivalent fits to 
the measured 7r° and tt^ [43| spectra indepen- 

dently for various centrality selections. Figure [TT] shows 
the comparison of data and parametrization on absolute 
and relative scales for minimum bias Au-|-Au collisions 
as an example. 

It is obvious that the neutral pion spectrum from 2004 
is systematically different from the 2002 result. The dif- 
ference between the parametrization and both sets of 
neutral pion spectra gives rise to an additional system- 
atic uncertainty in the resulting electron cocktail. The 
difference between the parametrizations is assumed to 
be one standard deviation of error due to this system- 
atic difference. This systematic uncertainty reaches its 
maximum value for electrons in the px range between 3 
and 5 GeV/c and is significantly less for lower as well as 
higher pT- 

Given that pion decays are the most important cock- 
tail ingredient at low and intermediate px it is obvious 
that the cocktail systematic uncertainty is largely domi- 
nated by the uncertainty in the pion spectra as well. To 
evaluate this uncertainty the full cocktail calculation is 
repeated with the pion cross section moved up or down 
by one standard deviation in the systematic uncertainty, 
propagating the uncertainty in the pion spectra to the 
electron cocktail. With a systematic uncertainty of 10 % 
almost independent of pt , some of which originates from 
the difference between the Run 2 and Run 4 7r° mea- 
surements, the pion input represents the largest contrib- 
utor to the electron cocktail uncertainty, except for at the 
highest Pt where direct radiation becomes important. 

Other light mesons contributing to the electron cock- 
tail are the 77, p, w, ry', and <f) mesons via their Dalitz de- 
cays and/or the conversion of photons from their decays. 
However, only the 77 meson is of any practical importance 
here. 

For the cocktail calculation, the shape of the invari- 
ant Pt distributions and the relative normalizations to 
the TT^ are required as input parameters. The pt spectra 
are derived from the pion spectrum by mr scaling, i.e. 
the same modified Hagedorn parametrizations are used 

(Eq. [6]), but with pT replaced by \Jp't + "^^meson ~ '^^o- 
The resulting jy/Tr*^ ratios agree as a function of pt well 
within experimental uncertainties for pt > 2 GeV/c with 
corresponding PHENIX data for p+p and Au-|-Au colli- 
sions. 



TABLE VI: Ratios of mesons to neutral pions in p+p collisions 

meson to pion ratios 
77/71" = 0.48 ± 0.03 [50] 
p/tt" = 1.00 ± 0.30 [51] 
lj/tt° = 0.90 ± 0.06 [52] 
■n'/ir° = 0.25 ± 0.075 [51J 
(j)/-K" = 0.40 ±0.12 [51] 



Since the chosen approach of rriT scaling ensures that 
at high Pt the spectral shapes of all meson distributions 
are the same, the normalization of the meson spectra 
relative to the pion spectrum can be given by the ratios 
of mesons to pions at high pT (5 GeV/ c is used here) . The 
values used for p+p collisions are shown in Table I VII 

For Au-j-Au collisions the same central values are used, 
but the uncertainties of the ratio rj/ir'^ and lu/tt^ are, con- 
servatively, increased to 0.10 and 0.27, respectively, since 
precision measurements are available for the production 
of these mesons in p+p collisions but not in Au-f Au col- 
lisions at all centralities. 

The contribution from the K^s decay and the semi- 
leptonic decay of Kg can only be determined via a full 
GEANT simulation, taking into account the exact elec- 
tron identification cuts. The electron cocktail includes 
parametrizations based on such GEANT simulations us- 
ing the charged kaon spectra measured in p+p and 
Au-f Au collisions as input [i^ . Systematic uncertainties 
are fully propagated from the kaon spectra to the spectra 
of reconstructed electrons from Kg^ decays. In any case, 
the contribution from kaon decays is only relevant (i.e. 
larger than 10 %) for electrons with pT < I GeV/c. 

Contributions to the electron cocktail from direct radi- 
ation are twofold. First, real photons produced in initial 
hard scattering processes, i.e. so-called direct photons, 
convert in material in the PHENIX aperture exactly as 
photons from light neutral meson decays. Second, every 
source of real photons also presents a source of virtual 
photons. Consequently, direct real photon production 
is accompanied by direct virtual photon production, i.e. 
the emission of e~'"e~ pairs. In the case of the neutral pion 
these two sources are the the 77 decay of the 7r° and the 
corresponding Dalitz decays, which are also called inter- 
nal conversions. The measured real direct photon spectra 
are parametrized and the conversion electron spectra of 
these are added to the electron cocktail. 

Figure [12] shows comparisons of the measured direct 
photon spectrum with the cocktail parametrization for 
(upper) p+p [ssj and (lower) minimum bias Au-|-Au col- 
lisions |54| . In accordance with direct photon measure- 
ments in Au-j-Au collisions, the direct photon yield is 
assumed to scale with the number of binary collisions as 
a function of the centrality for pT below 12 GeV/c in 
the Au-j-Au electron cocktail. Note that although the 
parametrization used underestimates the central value of 
the direct photon measurement below 5 GeV/c, the con- 
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FIG. 12: (Color online) Measured direct photon spectrum 
(red squares) compared with the cocktail parametrization 
(histogram) for (upper) p+p and (lower) Au+Au collisions. 



tribution from direct photons in this pr range is small 
enough that effects of this slight underestimation can be 
neglected. 

The ratio of virtual direct photons to real direct pho- 
tons depends on px because the phase space for dielec- 
tron emission increases with increasing pT [ssj . The very 
same effect is seen in the Dalitz decays of light neutral 
mesons. For instance, the Dalitz decay branching ratio 
relative to the two photon decay branching ratio is larger 
for the r] meson than for the 7r°. Consequently, the ratio 
of virtual and real direct photon emission increases with 

PT- 

Figure [13] shows the resulting electron cocktails for 
p+p and Au+Au collisions at 200 GeV. Systematic er- 
rors are estimated for all cocktail ingredients, propagated 
to the corresponding electron spectra, and then added in 
quadrature to determine the total cocktail systematic er- 
ror. The following systematic errors are assigned to the 
various inputs: 

• pion spectra: obtained via full cocktail calculations 
using pion spectra moved up (down) by the sys- 
tematic uncertainty of the pion spectra as input (al- 



most no pt dependence). This is the dominant sys- 
tematic uncertainty, except for the case of high pT 
electrons (above 5 GeV/c) in central Au-t-Au colli- 
sions, where direct radiation dominates the electron 
background and, at the same time, dominates the 
systematic uncertainty in the electron cocktail. 

• meson to pion ratios: the systematic uncertainties 
are listed in Table IVTl Since the contributions from 
all other mesons are much smaller than the contri- 
bution from 7] decay only the rj is of any practical 
relevance. This contribution is small compared to 
the uncertainty in the pion spectra and it depends 
only slightly on pT- 

• conversion material in the aperture: the contri- 
bution from photon conversions obviously depends 
on the material present in the aperture. A care- 
ful analysis of fully reconstructed dielectrons from 
photon conversions suggests that this uncertainty is 
not larger than 10 %. It is crucial to note that the 
contribution from photon conversion to the back- 
ground electron spectra is less than half of the 
contribution from direct Dalitz decays as a conse- 
quence of the careful minimization of the material 
budget in the PHENIX central arm aperture. For a 
reliable measurement of nonphotonic electrons this 
is an essential prerequisite. 

• i^es decay: this contribution is estimated via a full 
GEANT simulation. Given the limited statistics of 
this calculation a 50 % systematic error is assigned. 

• direct radiation: this contribution is directly prop- 
agated from the systematic error quoted for the di- 
rect photon measurement. It is relevant only at 
high Pt in central Au-)-Au collisions. 

In previous analyses, the single electron background 
cocktail has not included contributions originating from 
quarkonia decay (J/^A and T) or Drell-Yan processes. 
Each of these processes has a small total cross section 
relative to the electron cocktail background. However, 
with increasing pT these processes begin to contribute 
meaningfully and are included in the cocktail. 

PHENIX has measured the J/'0 pr spectrum from 
to 9 GeV/c hip+p collisions [53 • To determine the ap- 
propriate input to the Monte Carlo event generator used 
to calculate the components of the electron cocktail, two 
functional fits to the J /^j data are performed. The first 

uses the so-called Kaplan function po (^1 + {vr/pif'^ , 
and the second assumes the rriT scaling function used in 
[57l |. Since both functions may provide valid representa- 
tions of the true J/^ spectral shape, the average of these 
two functional fits to the data determines the central val- 
ues used in the cocktail. The upper (lower ) systematic 
uncertainties of the spectral shape are determined by nor- 
malizing each of the functional forms up (down) by 10% 
and using the largest (smallest) result for the systematic 
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FIG. 13: (Color online) Invariant differential cross sections and multiplicity of electrons for electron cocktails in p+p and 
Au+Au collisions at 200 GeV for the indicated centrality ranges. 



error. In practice, the lower bound of the spectral shape 
is determined by the smaller Kaplan function normalized 
down a further 10%, while the upper bound is set by the 
larger tot scaling function normalized up 10%. 

While measured J/4> Pt spectra exist at RHIC ener- 
gies, the T Pt spectra have not been measured. Since 
the overall production cross section for T is estimated to 
be roughly 1% of that of the J/ip, contributions from T 
decay contribute much less to the single electron cock- 
tail. A NLO calculation for T production in the color 
evaporation model provides the input to the Monte Carlo 
calculation of the central value estimates fssi] . Without a 



measured T px spectrum to characterize the systematic 
uncertainty to the Monte Carlo input, the same relative 
systematic uncertainties derived from the J/ip data are 
applied to the T cocktail estimate. Compared to the 
quoted uncertainties in the NLO calculation, the use of 
data-driven systematic uncertainties for the T provides 
a more conservative estimate. 

A leading-order (LO) Drell-Yan calculation [H^ of the 
single-inclusive lepton spectrum, p-l-p {e~^ +e~)/2+ 
X, is also included in the updated electron cocktail. The 
calculation is for < 0.5 and uses the CTEQ6M parton 
distributions [13]. The scale used is pT, no cut is placed 
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FIG. 14: (Color online) Electron background cocktail with 
quarkonium and Drell-Yan contributions for p + p collisions. 



FIG. 15: (Color online) Electron background cocktail with 
quarkonium and Drell-Yan contributions for 0-20% centrality 
Au+Au collisions. 



on the lepton pair mass, and a K-factor of 1.5 is applied. 

Figure [Til shows the cocktail for p+p collisions with the 
quarkonium and Drell-Yan contributions. The bottom 
section of the plot shows the ratio of the cocktail with 
the quarkonium and Drell-Yan to that without. 

PHENIX has measured the J/^ px spectrum out to 
5 GeV/c in Au-hAu collisions at 200 GeV [6l|. In order 
to extrapolate out to higher px, the measured J/^ spec- 
trum from p + p collisions is scaled by Raa and A^'coU • 
Above 5 GeV/c, two extreme scenarios are considered, 
and the difference in the two resulting estimates for the 
extrapolation is assumed as a systematic error. In the 
first scenario, the Raa is kept constant from its value at 
5 GeV/c. In the second scenario, the Raa is assumed 
to increase linearly from its value at 5 GeV/c up to a 
value of 1 at 10 GeV/c, above which it is assumed to be 
constant. Figure [15] shows the cocktail with quarkonium 
and Drell-Yan contributions for 0-20% centrality Au+Au 
collisions. 

The J/'^ spectrum in Au-)-Au collisions has not been 
measured separately for 0-10% and 10-20% centralities. 
In order to estimate the electron background in these 
centralities separately, we assume that the Raa of J/^ 
is the same in 0-10% and 10-20% centrality collisions. 
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FIG. 16: Invariant yields of inclusive electrons with (open 
circles) and without (filled circles) the photon converter for 
MB events. 

These yields can be expressed as the following rela- 
tions: 



J. Converter subtraction method 

The yields of photonic and nonphotonic electrons are 
obtained by measuring the difference between inclusive 
electron yields with and without a photon converter of 
precisely known thickness: a brass sheet of 1.680 % ra- 
diation length (Xq). Figure [TCI shows the corresponding 



^Conv-out ^ ArT'+7vNon-7^ (7) 

^Co„v-in ^ ij^iv; + (1 - e)7V^^°"-'^. (8) 

Here, TVf (7vf°"^^°"t) is the measured electron 

yield with (without) the converter. {N^°'^~'^) is the 
photonic (nonphotonic) electron yield, e represents a 
small loss of 7vNon-7 ^.j^g converter. This blocking 
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factor has been aheady evaluated in a previous measure- 
ment (e ~ 2.1%) [2^. The main issue in this calculation 
is to determine i?-^, i.e, to understand how much the pho- 
tonic electron yield is increased by the converter. The 
main source of photonic electrons is a mixture of mesons 
(tt*^, 1], 77', Lu, and (fi) decaying into real or virtual photons 
with their different pr slopes. 



TABLE VII: Radiation length (L) of material near the inter- 
action point. The conversion probability (pConv^ calculated 
for the case of electrons emitted at =1.0 GeV/c. 



Material 


L (Xo) 


^Conv 


Beam pipe (Be) 


0.288% 


0.201% 


Air (r < 30 cm) 


0.099% 


0.069% 


Total 


0.387% 


0.270% 


Converter (brass) 


1.680% 


1.226% 



To calculate R^, it is necessary to know exactly the 
material amounts near the interaction point. Table lylll 
shows a list of each material thickness in units of radi- 
ation length. To reduce the air contribution, a helium 
bag was installed in the space surrounded by the DC 
(0.3 < r < 1.0 m). The converter sheet was rolled just 
around the beam pipe in converter runs. The conversion 
probability (pConv^ Tab. lVIIl is calculated for the case 
of electrons emitted at Pj^ = 1.0 GeV/c. The equivalent 
conversion probability of a virtual photon in tt" Dalitz 
decay (pDaiitz-j jg o.598% R-y can be estimated with 
these values at p^ = 1.0 GeV/c. 



JjCon 

^Conv 



:>Dalitz 



(with converter) 



pDahtz (-ftrithout converter) 



2.41. (9) 



To obtain a more realistic value of R-y considering geo- 
metrical effects, GEANT-based Monte-Carlo simulations 
[46j for photon conversions were performed with and 
without the converter. R,y is determined for tt^ and rj sep- 
arately. We use the 7r° spectrum measured by PHENIX 
as the input for the tt^ simulation and assume uit scal- 
ing (pt \Jpt ^ ^ "^^oi normalized at high px to 

ry/TT^ = 0.48 ± 0.1) to obtain the input for the ry sim- 
ulation. Since the r/ mass is larger than the 7r° mass, 
the phase space of 77 Dalitz decay is slightly larger than 
that of tt" Dalitz decay. The relative branching ratio 
(Dalitz decay)/(two 7 decay) is 1.2% for 7r° and 1.5% 

for 77 [gS^I. This difference makes smaller than R^ . 
Contributions from other mesons which undergo Dalitz 
decay (77', p, w, and 0) are smaU (6% at pT = 3 GeV/c, 
and smaller at lower px)- The particle ratios at high pr 
{rj'/TT° = 0.25 ± 0.13, p/7T° = w/ttO = 1.0 ± 0.5, (p/n^ = 
0.4±0.2) are used to determine R^. The ^/7r° ratio used 
here is consistent with our ir^ and (f> measurement. The 
uncertainties in the particle ratios are included in the 
systematic uncertainties of i?-^. For this method, it is es- 
sential that the amount of material is accurately modeled 



in the simulation. We compared the yield of identified 
photon conversion pairs in the data and in the simula- 
tion, and concluded that the simulation reproduces Rj 
within ±2.7%. This uncertainty is included in the overall 
systematic uncertainty. 

In the upper plot in Fig. \17\ R-y is indicated as a solid 
curve, which is compared with the ratio of the inclusive 
electron yield with/ without the photon converter {Rcn)- 



RcN = 



N. 



Conv — in 



N, 



Conv— out 



R^, + {1- e) j?NP 
1 + -Rnp 



(10) 



Here, i?NP is the ratio of nonphotonic/photonic electron 
yields {N^°''^~''' /N^)- If there were no nonphotonic con- 
tributions (i?NP = 0), then we would have i?cN = ^7- 
The lower plot in Fig. [T7] shows that i?cN gradually de- 
creases with increasing pj,, while R^ slightly increases 
with p^. The difference between i?cN and R-y proves the 
existence of nonphotonic electrons. 

The converter method is applied for 0.3 < p^ < 1-6 
GeV/c, and the cocktail method is applied for pf^ > 1.6 
GeV/c. Boxes in the lower plot in Fig. [17] are systematic 
errors of i?NP from each method. i?NP increases with 
Px and is more than 1.0 in p?p > 1.6 GeV/c. This plot 
gives an important demonstration that the amount of 
conversion material is relatively small in PHENIX. 

To tune the normalization of the photonic electron 
spectra in the cocktail, the converter/cocktail ratio of 
photonic electrons is calculated for all centrality classes 
in Au+Au collisions. The ratio in MB is fitted by a con- 
stant (converter/cocktail — 1.182). The constant is ap- 
plied to correct the cocktail spectra of photonic electrons 
for all centrality classes in Au-|-Au collisions. The ra- 
tios of photonic electron yield from the converter method 
to that from the corrected cocktail method for all cen- 
trality classes in 2004 Au-|-Au collisions are shown in 
Fig. [iSl The curves in each plot indicate the systematic 
errors of the cocktail method. Figure [19] shows the con- 
verter/cocktail ratio (after scaling by 1.182) of photonic 
electrons in 2005 p+p collisions. 



K. Kaon decay background 

After the subtraction of photonic electrons, almost all 
of the remaining background is from kaon decays. Kaons 
contribute to the nonphotonic electron spectra via Kc3 
decay and are only substantial at low p^ {< 1 GeV/c). 
Since kaons have a long lifetime, the decay electrons do 
not originate from the primary collision vertex. There- 
fore, their contribution can only be determined reliably 
in a full detector simulation. This simulation was per- 
formed similarly to the ir^ simulation. It is done for 
charged kaons and K1 to estimate the final electron rate. 
The input spectra have a parametrization obtained by 
the PHENIX kaon measurements [49] . 

To compare kaon contributions with photonic and non- 
photonic electrons, the ratio of the electron yield from 
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FIG. 17: (Color online) (upper) Ratio of inclusive electrons 
with/without the converter {Rcn, points with systematic 
error boxes) and ratio of photonic electrons with/without 
the converter {R-,, solid line with a systematic error band), 
(lower) The ratio of nonphotonic yield to photonic electron 
yield {Rnp, filled circles) and the ratio of electron yield from 
kaon decays to photonic electron yield {Rkp, filled squares) 
as a function of for MB. Filled circles with statistical er- 
ror bars indicate -Rnp produced by both the converter and 
cocktail methods. 



is the product of the geometrical acceptance and recon- 
struction efficiency described in IIIIFI For the PH data 
set, Aerec also includes the PH trigger efficiency described 
in IIII Gl The cross sections from the MB and the PH 
data sets are consistent with each other in overlapping 
Pt regions. 

M. Systematic uncertainties in invariant spectra 

In the p + p analysis, systematic uncertainties are cat- 
egorized into (a) inclusive electron spectra, (b) cocktail 
subtraction, and (c) converter subtraction. Category (a) 
is common to both the converter and cocktail methods, 
and includes the uncertainties in luminosity (9.6%), geo- 
metrical acceptance (4%), elD efficiency (3%) (from com- 
paring to conversion electrons and simulation), and the 
PH trigger efficiency (3% at the plateau). Uncertainties 
in cocktail subtraction (category (b)) include the normal- 
ization (8%) and px dependent shape uncertainty (2% at 
Pt — 2 GeV/c, increasing to 6% at 9 GeV/c) which orig- 
inates from increased measurement errors on the back- 
ground sources. In the converter analysis (category (c)) 
the dominant uncertainties are in Rj (2.7%) and in the 
relative acceptance in the converter and the normal runs 
(1.0%). 

In the Au+Au analysis there is no trigger efficiency 
correction, but there is additional uncertainty from the 
correction due to multiplicity-dependent efficiency loss. 
A conservative systematic uncertainty of 4% was added, 
ascertained from the differences of the two columns in 
Table 13 

These uncertainties are propagated into the uncer- 
tainties in the heavy-flavor electron yields and added 
in quadrature. At low (high) pT the uncertainties are 
amplified (reduced) to large (small) uncertainties in the 
heavy-flavor electron signal by approximately a factor of 
1/{S/B), where S/B is the ratio of the nonphotonic elec- 
tron yield to the photonic electron yield. 



kaon decays over the photonic electron yield (Rkp) is 
shown as filled squares in the lower plot of Fig. [T71 The 
background from kaons becomes negligible above a p^ of 
1 GeV/c. 

L. Invariant spectra 

The invariant cross section for single electrons from 
heavy-flavor decay in p + p collisions is calculated using 
the following formula, 

^dVHF _ 1 1 ^HF 1 ^^^^ 

dp^ C 2ttpt ApxAy Aercc Cbias ' 

where C is the integrated luminosity; A^^^ is the electron 
yield from heavy-flavor decay after subtraction of pho- 
tonic and nonphotonic background contributions; Aercc 



N. Inclusive electron V2 

The value of the inclusive electron V2 was measured 
with the reaction plane method which can be written as 

^ =Ao{l + 2i;2cos(2(0 -*))}, (12) 

where Nq is a normalization constant, is the azimuthal 
angle of electrons, and ^ is the reaction plane angle. 
The reaction plane was determined from the multiplic- 
ity in each segment of the beam-beam counters. As the 
measurement of the reaction plane is sensitive to non- 
flow effects such as jets, resonance decays and HBT cor- 
relations, the reaction plane was measured in two well- 
separated rapidity intervals. Since each BBC is roughly 
three units of pseudorapidity away from the central arms, 
non-flow effects in the reaction plane measurement are 
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FIG. 18: (Color online) The ratio of photonic electrons measured by the converter method to the cocktail calculation in 
Au+Au collisions for the indicated centralities. The upper and lower curves show the systematic error of the cocktail. Error 
bars are statistical only. 
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FIG. 19: (Color online) The ratio of photonic electrons mea- 
sured by the converter method to the cocktail calculation in 
2005 p+p collisions. The data from the MB (PH) data set 
are shown below (above) 1.8 GeV/c. The upper and lower 
curves show the systematic error of the cocktail. Error bars 
are statistical only. 



higher order Fourier coefficients do not contribute much 
to the measured azimuthal electron distribution. Using 
the reaction plane measured with the BBCs, the inclusive 
electron V2 was calculated as V2 = (cos(2(0 — ^'moas)))- 
The true V2 with respect to the true reaction plane can 
be expressed in terms of the observed V2 with respect to 
the measured reaction plane as [H, [65j 



V2 



(C0s(2(^'r 



*true))) ■ 



(13) 



Two methods were used to extract (cos(2(^nica 
^'truo)))- The first is an analytical calculation [M. IgI 



(cos(2(*. 



^'truo))) 



(14) 



expected to be small. Figure [20] shows the azimuthal 
distribution of inclusive electrons for various centrality 
and pt ranges, fit to the form of Eqn. [T^l One can see 
from the agreement of the curves and data points that 



2^2 



Xexp(- 



4 n^oi + hi f)] 



where x is equal to v'^°'^^V2N, and N is the BBC multi- 
plicity. This method was verified to be well approximated 
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FIG. 20: (Color online) Azimuthal distributions relative to 
the reaction plane for various centrality and transverse mo- 
menta. The distributions have been normalized to their aver- 
age value and centered around 0. 
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FIG. 21: Reaction plane resolution as a function of centrality. 

by [el 

(cOs(2(5'moas - *truo))) ^ 

with ^'mcas is the measured reaction plane using only the 
north (south) BBC. 

Figure [21] shows the reaction plane resolution deter- 
mined from Equation [151 Less than 10% of background 
remains due to accidental RICH associations of hadrons. 
Contributions from such background sourced was sub- 
tracted as: 



2(cos(2(^P^_ - 



))> 
(15) 



dN 



dN": 



back 



dicjj-'i'R.P.) d{(j)-^R,p,) d{(j)-^B..P.) 



(16) 



where N^^^^j^ is the number of electrons identified by the 
RICH and N^^^^ is the number of the background parti- 
cles. The number of random associations are obtained by 



replacing the north and the south sides of the RICH in 
the software, as described previously. The V2 for a given 
centrality bin [a, 6] can be expressed as 



.bin 



6 „m°as(^) 



dN^ 
dC 



res(C) 



-dC 



I 

J a 



^^^^ 



(17) 



where res{C) is the reaction plane resolution for a given 
centrality C. Due to lack of statistics, Equation [T7l was 
approximated by a Riemann sum over 10% centrality 
bins. It was verified that V2 is not sensitive to changing 
the bins for the Riemann sum. This is because the inte- 
grand in the numerator of Equation [T7| is quite flat with 
centrality, as the resolution grows with both the mea- 
sured V2 and the particle multiplicity. Figure [521 shows 
the transverse momentum and centrality dependence of 
the inclusive electron V2- The minimum-bias V2 is shown 
in Fig. m 



O. Background V2 cocktail 

The azimuthal distribution of electrons {dN'^/d(t)) is 
given as the sum of the azimuthal distributions of 
photonic electrons {dN'' /dcf)) and nonphotonic electrons 
((iiV"°"-T/d(/)): 



dN7 dNT'^'i 



(18) 



Expanding Eq. [THjinto the Fourier expansion, the non- 
photonic electron V2 can be expressed as 



(1 + -RNp)f2, 



i?NP 



(19) 



where V2^ is the V2 of inclusive electrons, is the V2 
of the photonic electrons, and i?NP is the ratio of non- 
photonic to photonic electron yields which was defined in 
section UlIJI The photonic electron V2 can be estimated 
from the V2 of electrons from various photonic sources as 



RiVn 



(20) 



where Ri is the relative contribution of electron source i 
to the background and w| is electron V2 from the electron 
source. 

The dominant sources of photonic electrons are photon 
conversions and Dalitz decays from tt*^. In addition, we 
also took into account electrons from 77 and direct photon 
decays when calculating the photonic electron V2- The 
relative contributions of electrons from those sources are 
shown in Fig. 1241 The other sources were ignored when 
calculating the photonic electron V2 due to their small 
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FIG. 22: pt dependence of inclusive electron V2, for different centrality bins. A 5% systematic scale uncertainty is not shown. 



contribution. The measured V2 for tt and tt'^ was used 
as input for a simulation of tt'^ e^V2, and transverse 
kinetic energy scaling J^] was assumed for the 77, p, and 
w. The direct photon V2 was assumed to be zero. The 
decay electron V2 from ir^ and rj from Monte Carlo sim- 
ulations are shown in Fig. [25] and Fig. [26l The photonic 
electron V2 that was calculated from the results is shown 
in Fig. [27l The middle solid line is the mean value of 



the photonic electron V2 and the upper and lower dashed 
lines show one standard deviation of the systematic un- 
certainty. The uncertainty mainly comes from the uncer- 
tainty of the measured input V2 ■ The V2 from photonic 
sources for various centralities is shown in Fig. [271 
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FIG. 23: (Color online) Inclusive electron V2 for minimum 
bias events as a function pr measured in Au+Au collisions at 
y^s^ = 200 GeV. 



FIG. 25: (Color online) Electron V2 from ttq decay as a func- 
tion of pt, from a simulation. 




4.5 5 

(GeV/c) 



FIG. 24: (Color online) Relative contributions of electron 
from tt", rj and direct 7 to the background. 



P. Photonic V2 converter method 

The photonic electron V2 can be also determined by the 
converter method. Nonphotonic and photonic electron 
V2 can be separated by using the inclusive electron V2 
measured with (w^""""*") and without („™«''-°"*) the 
converter as 



non— 7 



conv — out 



conv — in 



Rnp{R^ — 1) 



(21) 



(22) 



Figure [55] shows the inclusive electron V2 with and with- 
out the converter. If the photonic electrons and non- 
photonic electrons have the same V2, the V2 measured 
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FIG. 26: (Color online) Electron V2 from decay as a function 
of pt, from a simulation. 



with and without the converter would be the same. Due 
to the small statistics of the converter runs, V2 measured 
with the converter has a large statistical uncertainty. The 
photonic electron V2 obtained by the converter method is 
shown as open squares in Fig. 1291 The result is consistent 
with the photonic electron V2 determined by the cocktail 
method within statistical uncertainty. 



Q. Heavy flavor V2 and systematic uncertainties 

The remaining background after subtracting photonic 
background is kaon decay as described previously. The 
electron V2 from kaon decays was also calculated by a 
Monte Carlo simulation assuming transverse kinetic en- 
ergy scaling and it was removed from the nonphotonic 
electron sample. The electron V2 from kaon decays was 
subtracted from the nonphotonic electron V2 as 



heavy 



Rknpv^. 



1 — Rknp 



(23) 
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FIG. 27: (Color online) Electron V2 from photonic sources as 
a function of pr and centrality. 
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FIG. 29: (Color online) Photonic electron 112 determined by 
two independent methods, the cocktail and the converter 
method. The lines are determined by the cocktail method and 
the squared points are obtained by the converter method. 



suring the inclusive electron V2 while changing the 
electron identification cuts from the standard cuts 
which were described in previous section. In this 
analysis we changed the cut parameters oi E/p, nO 
and /npeO. The relative change in V2 from vary- 
ing the cuts on E/p is about 2%, while for nO it 
is about 2% and for /npeQ it is 1%. The total 
systematic uncertainty assigned from electron iden- 
tification is 3 %. 
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FIG. 28: (Color online) Inclusive electron «2 with and without 
the converter installed. 



Here Rknp is the ratio of the yield of electrons from 
kaon decays and from all other nonphotonic sources 

{Rknp = iVf /A^^"™"''), and ^2^™^ ^^e V2 of electrons 
from heavy flavor decays. After subtracting kaon decays, 
the main source of nonphotonic electrons remaining is 
heavy flavor decays. 

Systematic uncertainties of the heavy flavor electron 
V2 are summarized below. 

• Reaction plane determination: The systematic un- 
certainty of the reaction plane determination was 
estimated by measuring the inclusive electron V2 
separately with the north BBC, the south BBC, 
and both the north and south BBCs combined. The 
maximum difference is about 5 % and we apply it 
as the uncertainty due to the reaction plane deter- 
mination. 

• Electron identification: The systematic uncertainty 
from electron identification was estimated by mea- 



• Background V2'. As described in the previous sec- 
tion, the uncertainties from the photonic electron 
V2 come from the tt*^ V2 and kaon V2, which have 
about 5 % uncertainties. We applied these values 
to the uncertainty of the photonic electron V2 ■ We 
also applied 5 % for the V2 from three-body kaon 
decays. 

• i?Np: The systematic uncertainty of i?NP comes 
from the uncertainties of inclusive electron spectra 
and the subtracted background spectra. The sys- 
tematic uncertainty of the inclusive electron spec- 
tra includes the uncertainties in the geometrical ac- 
ceptance, the reconstruction efficiency and the oc- 
cupancy correction, as described previously. 

The total systematic uncertainty was obtained by a 
quadratic sum of the above uncertainties. 

IV. RESULTS 

A. Heavy-flavor electron cross sections (p -|- p) 

Fig. [30] (a) shows the invariant differential cross sec- 
tion of electrons from heavy-fiavor decays. Fig. [30] (b) 
shows the yield divided by a FONLL calculation . All 
background has been subtracted. The data from the two 
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analysis methods (converter and cocktail) are combined: 
at low pt {pt < 1.6 GeV/c) the converter subtraction 
method is applied to the MB data set; at intermediate pr 
(1.6 < Pt < 2.6 GeV/c) the converter method is applied 
to the PH data set; and at highpT {pt > 2.6 GeV/c) the 
cocktail method is applied to the PH data set. 

Figure I5l1 shows the ratio of nonphotonic electrons (in- 
cluding charmonium, bottomonium, and Drell-Yan) to 
the photonic background. 
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FIG. 30: (Color online) (a) Invariant differential cross sec- 
tions of single electrons as a function of pT in p + p collisions 
at ^/s = 200 GeV. The error bars (bands) represent the sta- 
tistical (systematic) errors. The curves are the FONLL cal- 
culations [g^. (b) The ratio of FONLL/Data as a function 
of Pt- The upper (lower) curve shows the theoretical upper 
(lower) limit of the FONLL calculation. In both panels, a 
10% normalization uncertainty is not shown. 



B. Heavy flavor electron invariant yield (Au+Au) 

The differential invariant yield spectra as a function of 
p^ for the measured signal electrons from heavy flavor 
decays are calculated as in Eq. [211 Nf^*^ is the number of 
events in a centrality class i {i = 0-10 %, 10-20 %, 20- 
40 %, 40-60 %, and 60-92 %). e^'^'°'^{p^ ) is th e hadron 
contamination factor that is mentioned in llllEI e^^(pf-) 
is the acceptance correction in llll Fl is the embed- 

ding efficiency in IIII HI Ap^ is the p^ bin width. Ay is 
the rapidity range {\y\ < 0.5) where the input are dis- 
tributed at the first stage of single particle simulation for 
the acceptance calculation fsee IIII F]) . NJ^^ {p^,e~) and 
NJ^^ {pfp,e'^) are the resulting counts of signal electrons 
and positrons from heavy flavor decays by the converter 
method or cocktail method. 
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FIG. 31: (Color online) Ratio of nonphotonic electrons to 
photonic back- ground. Error bars are statistical errors and 
the error bands show the cocktail systematic errors. The solid, 
dashed, dot- dashed, and dotted curves are the remaining 
nonphotonic back- ground from Ke3 , p — > e^e~, oj — > e^e~ 
, and hadron contamination, respectively. 
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Figure [31] shows the invariant yield of inclusive and 
nonphotonic electrons in Au+Au collisions for various 
ranges in centrality. From top to bottom the spectra cor- 
respond to data in MB events, the five centrality classes, 
and the 2005 p+p data. The various curves are fits to the 
p+p data scaled by powers of ten for clarity (each factor 
is shown in the legend of Fig. l32|) . These spectra are pro- 
duced by the converter method for 0.3 < pT < 1-6 GeV/c 
and by the cocktail method for 1.6 < pt < 9.0 GeV/c. 
The boxes and bars are systematic and statistical errors 
for each data point, respectively. Most of the nonpho- 
tonic electron yield is from the decay of open heavy flavor 
mesons. The curves overlaid in the right panel of Fig. [351 
are the flt to the corresponding data from p+p collisions 
with the spectral shape taken from a FONLL calcula- 
tion [6^ and scaled by the nuclear thickness function 
TauAu for each centrality class. They are also scaled by 
powers of ten for clarity. Each scaled FONLL curve al- 
most agrees with the measured data in the low-p^ re- 
gion. On the other hand, we can see the disagreement 
in the high-p^ region between the measured data and 
the FONLL curve in the minimum bias, 0-10%, 10-20%, 
and 20-40% centrality classes. The data points are lower 
than each FONLL curve at high-pT- To quantify the sup- 
pression, the nuclear modiflcation factor is calculated in 
the next subsection. 

The insert box in the right panel of Fig. 13^ shows the 
ratio of signal (electrons from heavy-flavor decays) to 
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FIG. 32: (Color online) Invariant yields of (left) inclusive electrons and (right) open heavy flavor electrons for different Au+Au 
centrality classes, scaled by powers of ten for clarity. Error bars (boxes) depict statistical (systematic) uncertainties. The insert 
shows the ratio of these electrons to those from background sources in minimum bias events. The curves are scaled fits to the 
p + p spectrum. 



background as a function of for minimum bias events. 
This signal/background ratio (i?sB) is calculated as: 



R 



SB 



N. 



HF 



KV 



N2 



(25) 



Here, ]\[^°'^-'y is the measured nonphotonic electron 
yield, A''^' is the yield of photonic electron background, 
and N^'^ is the yield of electron backgrounds from kaons, 
vector mesons, quarkonia, and Drell-Yan. Filled circles 
with statistical error bars show R^^ produced by both 
the converter and cocktail methods. Boxes are system- 
atic errors of Rsb- The yields from the converter and 
cocktail methods are combined at pT = 1.6 GeV/c. The 
signal to background increases rapidly with pT, reaching 
one for px ~ 1.5 GeV/c, reflecting the small amount of 
conversion material in the detector acceptance. 



C. -RauAu (pt) 

As seen in Fig. [32l the Au+Au spectra agree well with 
the p+p reference at low p^ for all centrality classes, but 
a large suppression with respect to p+p develops towards 
high p^ in central collisions. To quantify the suppression 



effect, the nuclear modification factors, RAuAuiPr) is cal- 
culated: 



-RauAu(Pt) 



(A^coii) X dN^p/dp^j 



AuAu 



X da^p/dp^ 



(26) 
(27) 



where dN^^j^^/dp^ is the differential invariant yield in 
Au-I-Au collisions and dNpp/dp^ {da^p/dp^) is the differ- 
ential invariant yield (cross section) in p+p collisions at a 
given bin. Forpf, < 1.6 GeV/c, dNpp/dpT (dapp/dp^) 
is taken from the 2005 p + p data. At higher pT, the fits 
shown in Fig. [32] are used to remove statistical fluctua- 
tions, and the statistical error is moved to a systematic 
error in the shape of -RauAu- 7auAu(&) is the nuclear 
thickness function for A-l-A at an impact parameter (6) 
and A^coii(&) = TauAu(&) x cr^pp- Here, tr™ is the inelas- 
tic scattering cross section of p+p collisions (42 [mb]). 
^AuAu(Pt) for each centrality class is shown in Fig. 1331 
Systematic errors consist of contributions from TauAu (a 
horizontal band around RauAu{Pt) = l-0)i and fromp+p 
and Au-|-Au data (vertical bands for each data point). 
The statistical errors are shown as vertical bars for each 
data point. Data points and errors for each centrality 
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class are listed in the Appendix. 

If no nuclear modification exists and binary scaling is 
correct, J?auAu should be unity. In all centrality classes, 
i?AuAu is consistent within errors with unity for pT < 2 
GeV/c. However, we can see very clear suppression for 
0-10%, 10-20%, 20-40%, and minimum bias events in 
the high-pT region. This observed strong suppression in 
the high-pT region indicates that the medium created by 
Au+Au collisions is so dense that not only light quarks, 
but also heavy quarks lose their energy, or merely that 
heavy quarks in Au+Au collisions are not produced as 
compared with scaled p+p collisions. The total inte- 
grated yield in Au+Au collisions is compared with the 
binary-scaled yield in p+p collisions in the next subsec- 
tion. 



1. The invariant cross section is extraploated to pt=0 
to obtain the integrated cross section da'^/dy at 
y=0 over all pr- 

2. The cross section is convertered to a charm cross 
section at y = {dacc/dyly^o) using a total electron 
branching ratio of charm (BR(c — > e)) as 



1 



1 da'' 



dac 



dy BR{c e) Ce/o dy 



(29) 



D. Integrated RauAu{N^ 



part } 



We calculate the p|.-integrated i?AuAu as a function of 



part 



as: 



-^AuAuC^pa 



NluAuiPT) 



(TauAu) X a^JpT) (^coll) X N!^Jpt) 



where N^^j^^{pt) is the total electron yield above a trans- 
verse momentum of pr- Figure l34l shows ^AuAu(^part) 
for electrons from heavy-flavor decays for six different 
integrated p^ ranges as a function of the number of par- 
ticipant nucleons, iVpart- When the lower limit of integra- 
tion is reduced to p^ = 0.3 GeV/c, which includes more 
than half of the heavy-flavor decay electrons predicted 
by the FONLL calculation in p+p collisions, i?AuAu is 
close to unity for all A'^pait- This behavior suggests that 
the total yield of electrons from heavy flavor decays in 
Au+Au collisions is the same as the binary-scaled yield 
in p+p collisions. The observed strong suppression phe- 
nomenon can be interpreted as the result of softened 
heavy quarks due to energy loss in the created medium. 
For p^ > 0.6, 0.8, 2.0, 3.0, and 4.0 GeV/c, i?AuAu 
decreases systematically with centrality. The behavior 
seems to be similar to the observed suppression for tt" 
and rj mesons [^,[5^. However, quantitative comparison 
with the suppression of the light mesons requires an un- 
derstanding of the decay kinematics of open charm and 
bottom mesons. 



E. Total charm cross section 

We determine the total charm cross section in p + p 
from the heavy-flavor electron invariant cross sections 
shown in Fig. [SO] The heavy-flavor electron cross sec- 
tion at low pt is domianted by charm, while the bottom 
contribution is small, order of 1%. 

We follow the following steps to determine the total 
charm cross section. 



Here BR{c e) is the total charm to electron 
branching ratio, and Ce/p — 0.935 is a kinemat- 
ical correction factor to account for the difference 
of rapidity distribution of electrons and D mesons. 

3. Finally, the cross section is extrapolated to the en- 
tire rapidity range to obtain the total charm pro- 
duction cross section CTcc- 

We use two methods for the step [TJ In the first 
method, the total charm cross section in p + p collisions 
is derived by numerically integrating the heavy-flavor 
electron cross section for pT > 0.4 GeV/c: dae{pT > 
OA)/dy = 5.79 + 0.59 + 1.64 ^b. Here the system- 
atic uncertainty is obtained by integrating the upper and 
the lower systematic uncertainty limits of the differential 
cross section, since the systematic uncertainties are es- 
sentially coherent. We don't include the lowest pT point, 
at Pt = 0.35 GeV/c since its large individual uncer- 
tainty would make the uncertainty in the numerical in- 
tegral larger. The cross section is then extrapolated to 
Pt — using the spectrum shape predicted by FONLL: 
dae{pT > 0)/dy = 10.6 + 1.1 + 3.2 ^b. Here we have 
assigned 10% (a quarter of estimated cross section for 
Pt < 0.4 GeV/c) to the systematic uncertainty of the 
extrapolation, and have subtracted contributions from 
b e and 6 — c — e cascade decays (0.1/ib). 

In the second method for the step [TJ we use PYTHIA 
to calculate the spectrum shape from charm decay elec- 
trons. Here we use PYTHIA as an event generator for 
_D-meson decays. Since electrons are from D and other 
charm hadron decays, its momentum distribution in low 
Pt part and the high pT part are correlated. If one can 
adjust the pt distributions of charm hadrons so that the 
electron spectrum for just above the un-measured region, 
Pt > 0.3 GeV/c, is reproduced, the pt distribution in the 
lower Pt region should be well described by the simula- 
tion. 

Figure 135] shows comparisons between the PYTHIA 
calculations and the data. The five curves in are from 
PYTHIA calculations, and the red points are from the 
data. The PYTHIA curves are fitted to the data for the 
fit range 0.3 < pt < 1.6 GeV/c, with their normalization 
to be the sole fitting parameter. The shapes of PYTHIA 
curves are fixed. The PYTHIA curves underestimate the 
data for high pT {pt > 1-6 GeV/c). However, in the 
Pt ranges that is used for the fit, the PYTHIA curves 
can reasonably describe the data. The five curves cor- 
responds to 5 different choices of < fc^ > ranging from 
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FIG. 33: (Color online) Open heavy flavor electron -RauAu for the indicated centralities. 



1.5 GeV/c to 3.5 GeV/c shown in the legend. The five 
curves are very similar in the low px region below the 
data. This is because the electron spectra in the region 
is mainly determined by the decay kinematics. Since the 
sohd black curve (< kr >= 2.5 GeV/c) describes the 
data well, we use it as the central value. We note that 
similar values of < /ct > are found in our two hadron 
correlation analyses in p + p [gtI ] . 



Using the second method, we obtained da^ / dy{pT > 
0) = 10.9 ± 0.6 ± 2.8/ib. Here we have subtracted beauty 
contribution and added the systematic uncertainty of 
±1.4yLtb in quadrature as the systematic uncertainty of 
this extrapolation procedure. This systematic uncer- 
tainty is estimated as 50% of the difference of the cross 
section with < kx >= 3.5 GeV/c (minimum cross sec- 
tion) and with < kx 1.5 GeV /c (maximum cross sec- 
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FIG. 34: (Color online) Nuclear modification factors -RauAu for open heavy flavor electrons vs centrality, integrated above the 
indicated Pt ranges. 



tion). Since the second method gives a value of da'^/dy 
very close to that of the first method, we use the value 
obtained from the first method in the charm cross section 
calculation. 

We determine the charm production cross section, 
daac/dy = 119 ± 12 ± 38 /xb using Eq. (gH). The branch- 
ing ratio, BR{c — > e) = 9.5 ± 1.0%, is calculated us- 



ing the following charmed hadron ratios: D^/D'^ = 
0.45±0.10, Ds/D° = 0.25±0.10, and Ac/D^ = 0.1±0.05. 
In the step |3l we calculated the rapidity distribution of 
charm quark using HVQMNR [H] with CTEQ5M [H 
PDF. The total charm cross section is determined to be 
a^c = 551 ± 57*^**1' ± 195''''^ ^b. Here we have assigned 
a systematic uncertainty of 15% to the extrapolation to 
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TABLE VIII: Charm cross section per N + N collision in centrality bins (0-10%, 10-20%, 20-40%,40-60%, 60-92% and Min 
Bias(MB)) in Au+Au and p + p. The second column showS the nuclear overlap integral Taa of the centrality. 



Centrality 


Taa (mb-i) 


%-\pT>o-i (^b) 


^ (Mb) 


dace 
dy 


(Mb) 


dec (/ib) 


0-10% 


22.8 ± 1.6 


6.50 ±0.12 ±0.92 


11.9 ±0.2 ±2.0 


134 ± 


2 ±27 


620 ± 11 ± 156 


10-20% 


14.4 ± 1.0 


6.30 ±0.13 ±0.96 


11.6 ±0.2 ±2.1 


130 ± 


3 ±27 


600 ± 12 ± 155 


20-40% 


7.07 ±0.58 


6.15 ±0.13 ±1.07 


11.3 ±0.2 ±2.2 


127 ± 


3 ±29 


586 ± 12 ± 159 


40-60% 


2.16 ±0.26 


6.64 ± 0.22 ± 1.28 


12.2 ±0.4 ±2.6 


137 ± 


5 ±33 


633 ± 21 ± 180 


60-92% 


0.35 ±0.10 


5.98 ± 0.42 ± 1.68 


11.0 ±0.8 ±2.4 


123 ± 


9 ±39 


570 ± 40 ± 201 


0-92% (MB) 


6.14 ±0.50 


5.96 ±0.08 ±0.96 


10.9 ±0.2 ±2.1 


123 ± 


2 ±27 


568 ± 8 ±150 


p + p 


NA 


5.79 ± 0.59 ± 1.64 


10.7± 1.1 ± 3.2 


119 ±12 ±38 


551 ±57 ±195 
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FIG. 35: (Color online) Comparison of electron pr distribuion 
calculated by PYTHIA and the data (filled circles). The five 
curves are electron pr distribuions from charm decay calcu- 
lated using PYTHIA with different < kr > values shown 
in the legend. All five curves are fitted to the data for 
0.3 <pT < 1.6 GeV/c. 



the full rapidity. 

The charm cross section obtained here is slightly 
smaller than that pubhshed in ^ (567 ± 57"*''* ± 
193''^'' /ib) but agrees well within the systematic uncer- 
tainties. The difference comes from the fact that we 
now subtract the contributions from and J/ip from 
the electron spectrum. The cross section is compati- 
ble with our previous measurement in the 2002 run (2^ 
(920 ± 150''*''* ± 540">"' fih) and agrees weU with the value 
derived from the yield of the dielectron continuum [st'I 
(544 ± 39"*''* ± 142"y" ± 200™°'^'=' fih). 

The FONLL cross section (256^^™ fib) is compatible 
with the data within its uncertainty. Although the data 
extends to high px where the bottom contribution is ex- 
pected to be dominant, the present analysis does not sep- 
arate charm and bottom contributions. The bottom cross 
section predicted by FONLL is 1.87t;67 A*b, and the up- 
per FONLL curve is consistent with the data. The bot- 
tom cross section by FONLL is also consistent with the 
bottom cross section in p + p we deduced from electron- 
hadron charge correlation [70| and dielectorn continuum 



mass distribution |57| . 

Total charm yields in Au+Au in various centrality bins 
are evaluated in a similar way from the heavy-flavor elec- 
tron spectra. For each centrality bin, the heavy-flavor 
electron spectrum is integrated for px > 0.4 GeV/c to ob- 
tain dN'^- /dy\p^yo,4- The electron cross section per N+N 



collision, da'^ / dy\p^yoA 



da'' 



is then obtained as 



1 dN" 



^lp.>0.4-^^|p.>0.4. 



(30) 



The cross section is then extrapolated to = and 
the contribution from b decays and b cascade decays are 
then subtracted to obtain the electron cross section from 
charm decay per N + N collision da^^/dy. Here we as- 
sume that Raa — 1-0 for pT-integrated b production. 
The cross section dacc/dy is obtained using Eg. ipS]) . Fi- 
nally, the cross section is extrapolated to the full rapid- 
ity range using the rapidtiy distribution of charm from 
HVQMNR to obtain total charm cross section acs per 
N + N collision. We use the same correction factors 
BR{c — > e), C^j-Q, and rapidity distribution as in p + p 
collisions. Table IVIIll summarizes the charm cross section 
per N±N collisions in each centrality bins of Au+Au col- 
lisions. The charm cross sections per + TV in Au+Au 
thus obtained agree well with the p+p cross section for all 
centrality bins. The data are consistent with the results 
from our previous measurement in Au+Au collisions in 
the 2002 run [26] (622 ± 57"*"* ± 160"^" \xh per NN). 

Figure [551 shows c?cr^g/dy|py>o.4 per N + N colhsion 
in Au+Au and p + p as a function of A^coU- The cross 
section is more than half (~ 54%) of the total electron 
cross section from charm decay da^^/dy. Fig. [33] shows 
that Raa is consistent with unity for pT < 1.4 GeV/c. 
Thus the shape of the electron spectrum in Au+Au and 
p + p is almost independent of Acoii in the low px region. 
This indicates that the fraction of electron cross section 
below pt < 0.4 GeV/c is almost independent of A'coii- 
In the following, we evaluate the Acoii dependence of the 
total charm cross section assuming that it is the same as 
the iVcoii dependence of dcr^g/d?/|py>o.4. 

To quantitatively evaluate the A^coU dependence of to- 
tal charm cross section, we fit the data with a power-law 
function A ■ iV^jj. We perform the fit with and without 
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TABLE IX: The results of power-law fits A ■ N^^^ shown in 
Fig. 1361 In the fit, only statistical error is considered. The un- 
certainty of A and the first error of a are the statistical erorrs 
obtained from the fit. The second error of a is systematical 
uncertainty. 





A{^^h) 


a 


xVndf 


Au-l-Au 


6.00 ±0.45 


1.0092 ±0.0120 ± 0.0506 


6.4/3 


Au-I-Au -1- pp 


5.98 ±0.34 


1.0097 ± 0.0094 ± 0.0403 


6.4/4 



the p + p data. The fitted curves are shown in Figure I36[ 
and the results of the fits are summarized in Table IIXI 

Most of the sources of systematic uncertainties in the 
charm cross section are independent of A'coii and they 
do not affect the value of a. For the systematic uncer- 
tainties in a, we consider the following three remainig 
sources: (i) Multiplicity dependence of electron recon- 
struction efficiency (ii) systematic uncertainties in Rj 
(iii) systematic uncertainties in Taa- The systematic un- 
certainty from these sources are propagated to a and 
added in quadrature to obtain the total systematic un- 
certainty, which is shown as the second error of a in ta- 
ble |lXl The value of a for Au±Au and p + p combined is 
a = 1.0097±0.0094"*''*±0.0403''y", consistent with unity, 
indicating that total charm cross section scales with A^coii- 
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FIG. 37: (Color online) Nonphotonic electron V2 as a function 
of pt in Au±Au collisions at ^/snn = 200 GeV. 

Therefore the nonphotonic electron V2 that was obtained 
by subtracting photonic electrons and kaon decays from 
inclusive electrons should be dominated by heavy fia- 
vor electrons, which refiects the azimuthal anisotropy 
of heavy quarks. In this momentum region, electrons 
mainly decay from D mesons. Therefore the non-zero 
nonphotonic electron V2 indicates that the D meson also 
has a non-zero V2- 

Figure [38] show the v^^ for various centrality ranges. 
The Pt binning used is coarser than that shown in the 
minimum-bias due to lack of statistics. The 
for 60-93% collisions is not shown as the pion V2 was 
not measured by PHENIX separately for that centrality 
range. At high px, the inclusive electron V2 in the 60- 
93% centraHty range was measured to be negative by 
between 2 and 3 standard deviations in the statistical 
error, and this accounts for the difference between the 
minimum-bias and the 0-60% . 

V. DISCUSSION 



FIG. 36: (Color online) A'coii dependence of charm cross sec- 
tion per N + N collisions in Au+Au and p + p collisions. The 
error bar show the statistical errors only. The curves in the 
figure is the power-law fit A ■ N^^n. 



F. Heavy flavor electron V2{pt) 

The nonphotonic electron V2 for minimum-bias events 
is shown in Fig. [37l The vertical lines are the statistical 
errors and the la systematic uncertainties are shown as 
bands. The nonphotonic electron V2 increases with pT 
below 2.0 GeV/c then it saturates or perhaps decreases. 
This trend is similar to that of the meson V2 (tt and K). 
The main source of the nonphotonic electrons is semilep- 
tonic decays of heavy fiavor mesons (charm and bottom) . 



A. Theory comparison (i?AuAu(PT)) 

Before the collection and analysis of RHIC data, it 
was generally expected that the dominant mechanism for 
the suppression of heavy quarks in the medium would be 
gluon radiation, and thus the heavy quarks would not 
be suppressed as much as the light quarks due to the 
dead-cone effect. Predictions from the BDMPS model 
and DGLV model for radiative energy loss can be seen 
in Fig. |39] j7l|, [t^] ■ These models do not account for the 
suppression of the heavy fiavor electrons, but they are 
sensitive to the relative admixture of D and B mesons 
produced in the collisions as D mesons are predicted to 
be suppressed due to gluon radiation more than are B 
mesons. However, recent measurements by PHENIX [t^ 
suggest that the admixture of D and B is comparable to 
what was used in the BDMPS and DGLV calculations. 
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FIG. 38: (Color online) v^^ for the indicated centralities. 



Mustafa [73[ found that radiative and elastic scatter- 
ing energy loss for heavy quarks are comparable over a 
very wide kinematic range accessible at RHIC. Contrary 
to what was previously thought, coUisional energy loss 
should be taken into account in the calculation of sup- 
pression of heavy flavor mesons in Au+Au collisions. Fig- 
ure [32] shows the DGLV prediction for suppression when 
coUisional energy loss is taken into account in addition 
to radiative energy loss Q. 

In 7J|, it was shown that a Langevin-based heavy 
quark transport model can qualitatively explain the large 
suppression (and azimuthal anisotropy) of electrons from 
heavy flavor decays in Au-|-Au collisions. The model 
places a heavy quark into a thermal medium, and as- 
sumes that the interaction of the heavy quark with the 
medium can be described by uncorrelated momentum 
kicks. Contrary to the models described above, the inter- 
action in the Langevin model is given exclusively by elas- 
tic collisions, which is a good approximation for quarks 
which are not ultra-relativistic in the center of mass 
frame of the collision. The parameter which is tuned in 
this model is the heavy quark diffusion coefficient. While 
the above Langevin model fails to simultaneously de- 
scribe Raa and V2 for a single value of the diffusion coef- 



ficient, another Langevin-based model [75|,|76[ is in good 
agreement with both the suppression and anisotropy. In 
this model, the elastic scattering is mediated by reso- 
nance excitation of D and B- like states in the medium. 
The theoretical evidence for the existence of such reso- 
nance states comes from lattice computations. Figure HOI 
shows the calculations from these two models. 

Recently, it has been suggested that coUisional disso- 
ciation of heavy quarkonia in the quark-gluon plasma 
fr^ may be a possible explanation for suppression of 
J/ip production in heavy ion collisions. A. Adil and 
I. Vitev investigated the pQCD dynamics of open charm 
and bottom production and, in the framework of the 
GLV approach extended to composite qq systems, de- 
rived the medium induced dissociation probability for 
D and B mesons traveling through dense nuclear mat- 
ter (tI]. They showed that the effective energy loss, 
which arises from the sequential fragmentation and dis- 
sociation of heavy quarks and mesons, is sensitive to the 
interplay between the formation times of the hadrons and 
QGP and the detailed expansion dynamics of hot nuclear 
matter. Figure |4T] shows their result as the band. 

Most of the models calculate the nonphotonic electron 
production assuming the same chemical composition of 
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FIG. 39: (Color online) -RauAu in 0-10% centrality class com- 
pared with energy loss models. The thick dashed curve is the 
BDMPS [tJ] calculation for electrons from D and B decays. 
The bands are DGLV ;73i calculations for electrons from D 
and B decays. The lower band contains coUisional energy loss 
as well as radiative energy loss. The thin dashed curves are 
DGLV calculations for electrons from D decays only. 
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FIG. 40: (Color online) Comparison of Langevin-based mod- 
els from [74|-[7g] to the heavy flavor electron -RauAu for 0-10% 
centrality and V2 for minimum-bias collisions. 

charm and beauty hadrons in p-l-p and Au-|-Au collisions. 
As it has been observed for the light hadrons [7S||, one 
could expect a modification of the charm hadron chem- 
ical composition in the most central Au-t-Au collisions. 
In particular, an enhancement of the Ac production has 
been predicted A Ac enhancement leads naturally 

to a nonphotonic electron Raa smaller than one due to 
a smaller semileptonic decay branching ratio of charm 
baryons compared to charm mesons and also due to a 
softer spectrum of the electrons from the charm baryon 



FIG. 41: (Color online) -RauAu in 0-10% centrality class com- 
pared with a coUisional dissociation model \7S\ (band) in 
Au+Au collisions. 



decay. A Raa of about 0.65 for electrons from charm 
hadron decay is predicted when a charm baryon to charm 
meson ratio in central Au-|-Au collisions close to one is 
assumed (sTj . 

Gossiaux and Aichelin [82| - [83 | calculated the Raa from 
coUisional energy loss in pQCD using a running coupling 
constant and replacing the Debye mass with a hard ther- 
mal loop calculation. The model finds a value close to 
the experimental Raa for all centralities, while leaving 
room for a possible radiative contribution as well. Fig- 
ure shows the Raa as a function of pt and centrality 
from this model. 

Kharzeev has predicted a universal bound on the en- 
ergy of a parton escaping strongly coupled matter (ssj in 
J\f — A SUSY Yang-Mills theory, under some assumptions 
about the evolution of gauge fields in heavy- ion collisions. 
In the model, the Raa at high px is given by a constant 
times A'part /NcoW with the constant universal for all par- 
ticle species. 



B. Theory comparison (112) 

The azimuthal anisotropy of heavy quarks gives infor- 
mation about the density of the colliding system. The 
collectivity of heavy quarks during the early stage of the 
collision is expected to be developed via the strong in- 
teraction inside the quark gluon plasma, which will be 
seen in the relatively low pT region. On the other hand, 
the suppression and/or modification of the highpT yield, 
which is caused by the interaction of the heavy quarks 
with the overlapped almond shape geometry and path 
length of the heavy quark inside the plasma, could also 
create an azimuthal anisotropy at higher px- Therefore, 
when considering the pr dependence of the anisotropy 
parameter V2, it is important to understand these dif- 
ferent effects as well as the relative fraction of bottom 
quarks compared to charm quarks. Several model pre- 
dictions are compared to the experimental data in this 
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FIG. 42: (Color online) -RauAu vs centrality and pr from Gossiaux and Aichelin [82148 



section. 

Figure |43] shows a comparison of two different model 
calculations to experimental data for the nonphotonic 
electron V2 from minimum-bias collisions. The models 
shown in this figure were also discussed in the previous 
section as they are also used to calculate Raa- In the 
model of Moore and Teaney [t^ , a Langevin-based heavy 
quark transport calculation is performed in a perturba- 



tive quark gluon plasma with hydrodynamic simulation. 
Two different calculations with two different assumptions 
for the value of the diffusion coefficient are shown as a 
solid curve for D = 3/(27rT), which agrees with experi- 
mental data, and as a dashed curve for D = 12/(27rT). 
In the model of Armesto et al. [7l| , the effect of heavy 
quark energy loss for pT > 2 GeV/c is considered without 
any collective effect in the lower pT region, which gives a 
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lower limit on V2 at high pT given by the energy loss and 
the geometry of the almond shape, with two values of the 
transport coefficient q for the green solid curve q — 14 
GeV^/fm and for the green dashed curve q — 4 GeV^/fm. 



• e- v""^ 

Moore SiTeaney 3/(2n:T) 

Moore STeaney 12/(27iT) 

Armesto et al. 14 GeV^/fm 
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FIG. 43: (Color online) Comparison of experimental V2 data 
with the models of Moore and Teaney [tJ] and Armesto et al. 

El. 



Figure |31] shows a set of comparisons of three model 
calculations with the experimental data. The model of 
Greco et al. [S^ is based on the quark coalescence model 
with light and charm quarks V2, where the light quark V2 
is estimated from the measured meson V2 and the same 
amount of radial flow is assumed for all quark species 
with the same saturating V2 for the charm quark as for 
the light quark. The calculated electron V2 from the D 
meson is labeled as "Greco et al. c flow" . On the other 
hand, the line labeled as "Greco et al. no c how" as- 
sumes zero V2 for the charm quark, and the predicted V2 
is given only by the light quark V2 via the quark coales- 
cence mechanism. The model of Zhang et al. [s^l is a hy- 
brid model starting with HIJING as an initial condition, 
followed by a parton cascade, and then a hadron cascade 
after a hadronization with a quark coalescence model. 
This is based on a multiphase transport model including 
rescattering of the charm quark with the other partons. 
The line labeled "Zhang et al. 10 mb" is the resulting 
electron V2 from D-meson decay with the charm quark 
parton-scattering cross section of Up = 10 nib, while the 
line labeled "Zhang et al. 3 mb" corresponds to ap — 3 
mb. The model of van Hess et al. [75,, J6.], discussed in 
the previous section as it is also used to calculate Raa, 
includes resonant interactions in a strongly interacting 
QGP and parton coalescence of both charm and bottom 
quarks. The band with shows the results of this model, 
while the result with only charm quarks is shown as the 



line labeled "van Hees et al. c only" , and the result with- 
out resonances is shown as the line labeled "van Hees et 
al. no reso" . Both the resonance and the coalescence 
are responsible for the large magnitude of the measured 
electron V2 in the low px region, while the reduction or 
flattening of V2 at high px requires the bottom quark 
contribution. 
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FIG. 44: (Color online) Theory comparison of experimental 
V2 data with the models of Greco et al. [s^, Zhang et al. [83], 
and van Hees et al. [tsI . [t^ . 

A pQCD-based elastic scattering model [82U83] with 
running coupling constant is used to calculate V2 as well 
as the Raa- The minimum-bias V2 from this model is 
shown in Fig. 1451 The f 2 is very sensitive to the interac- 
tion time, and a later freeze-out can produce a larger V2 




[GeV/c] 



FIG. 45: (Color online) Minimum-bias V2 cornpared with a 
model calculation from Gossiaux and Aichelin [82-[84ll. 



Various models discussed in this section indicate strong 
interaction of heavy quarks with the high energy density 
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matter created in heavy-ion collisions at RHIC and also 
support quark coalescence(recombination) as a dominant 
hadronization mechanism, which would mean that the 
heavy quarks are also a part of the quark gluon plasma. 
Therefore the sizable event anisotropy parameter V2 of 
electrons from heavy quarks measured at low pT should 
be related with the collectivity of the heavy quark. The 
reduction or flattening of V2 in the high px region would 
have to be caused by the change of the dominant mech- 
anism from the collectivity of the heavy quark to the 
energy loss of the heavy quark, especially when most of 
the heavy quarks are charm quarks. However the reduc- 
tion can also be qualitatively understood by an increase 
of the relative fraction of bottom quarks as a function of 
PT and a smaller magnitude of event anisotropy for the 
heavier quark. 

In order to determine the charm quark V2, the mea- 
sured nonphotonic electron V2 is compared to that ex- 
pected from coalescent production of D mesons. Chi 
square values are calculated for different combinations 
of light and charm quark V2 values. The pT range is 
restricted to below 2 GeV/c in order to be sensitive to 
only the charm quark, and not to the bottom quark. The 
quark coalescence is assumed to occur at similar veloci- 
ties, therefore it gives the momentum ratio of about 1/6 
for light/charm quarks, which is given by the ratio of the 
effective mass of the quarks. A common shape of V2 vs pt 
for both light and charm quarks is assumed, where this 
shape is determined by the light hadron V2- For given V2 
values of light and charm quarks, tt, K and proton 's 
and electron V2 from D mesons are calculated, then the 
total chi-square sum is calculated and plotted in Fig. 021 
as a function of the light quark V2 (horizontal axis) and 
charm quark V2 (vertical axis). Both the horizontal and 
vertical axes are normalized to the measured light quark 
W2, so that unity corresponds to the light quark V2- The 
chi-square minimum lies at about equal light quark and 
heavy quark V2 magnitudes, although the relative error 
in V2 is much larger for the charm quark than for the 
light quark. This indicates a common quark collectivity 
in the quark gluon plasma phase. 



C. Viscosity to entropy ratio 77/s 

The ratio of shear viscosity 77 to entropy density s is 
a key parameter that determines the damping rate in 
a relativistic system. At temperature T, characteristic 
damping times r are of order r]/sT fss!]. 

While the apparent success of ideal hydrodynamics in 
describing particle spectra and elliptic flow patterns at 
RHIC would imply a vanishing value of 77, straightfor- 
ward arguments based on the uncertainty principle sug- 
gest that the viscosity for any thermal system must be 
non-zero [1^. This observation was extended by Kovtun, 
Son and Starinets (KSS) [13, who demonstrated that 
conformal field theories with gravity duals have a ratio 
of viscosity 77 to entropy density s of l/47r (in natural 




u quark 

FIG. 46: (Color online) map as a function of the light 
quark V2 (horizontal axis) and charm quark V2 (vertical axis), 
each divided by the measured light quark V2 ■ 



units). KSS conjectured that this value is a bound for 
any relativistic thermal field theory, that is, rj/s > l/47r. 

The results presented in Section |IV] for RauAu{Pt) 
and V2{Px) show that heavy quarks lose energy in the 
medium, while acquiring a substantial component of the 
medium's collective flow. Both of these effects may be 
regarded as the damping of the initial non-equilibrium 
dynamics of the "external" heavy quark by the medium. 
The simultaneous description of these phenomena by the 
model of van Hees et al. provides a straightforward, al- 
beit indirect, method to infer the ratio of viscosity to en- 
tropy density. The resonance model employed in Ref. [75| 
leads to an estimate for the heavy quark spatial diffu- 
sion constant Dg ~ (4 — 6)/2ttT for temperatures T in 
the range 0.2 GeV < T < 0.4 GeV(see also Figure 23 
of Ref. [9ll|). Moore and Teaney [t^ perform a pertur- 
bative calculation of this quantity, and find that for a 
medium with three light flavors the ratio of Ds to the 
hydrodynamic diffusion constant r]/{e-\-p) for the bulk (e 
is the energy density and p the pressure of the medium) 
has a value of ^ 6 roughly independent of the coupling 
strength ~ mu/T, where mo is the Debye mass. They 
argue that the weak variation with coupling strength is to 
be expected in this ratio of transport coefficients, making 
it plausible that it remains near 6 in the strongly-coupled 
regime. In this case, and approximating the thermo- 
dynamic identity e + p = Ts + ^13^3 ~ Ts appropri- 
ate for the baryon-free central region, one readily finds 
rjl s ^ (1.33 — 2)/47r, that is, a value near the KSS bound 
and consistent with other estimates for the RHIC plasma 
based on flow ^] , fluctuations [11] , entropy produc- 
tion [il] and detailed hydrodynamic calculations [H, [o^] . 
It should be noted that these various estimates are based 
on observables of the "bulk" medium, in flavor channels 
dominated by u, c? and s quarks, while the result pre- 
sented here relies explicitly on the coupling of heavy fla- 
vor to the medium. The consistency of the derived value 
of rj/ s supports both the strong coupling of heavy flavor 
to the medium and the low value of ry/s for the RHIC 
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This article has detailed the measurement of the yield 
of single electrons from scmilcptonic decays of heavy fla- 
vor mesons at midrapidity in p + p and Au-|-Au collisions 
at -y/sNN = 200 GeV, as well as the azimuthal anisotropy 
parameter V2 of such electrons in Au-|-Au collisions. The 
unexpectedly large suppression of heavy flavor electrons 
in Au-I-Au collisions relative to those from p+p collisions 
and the large V2 of heavy flavor electrons has generated 
much theoretical work. In a system as complicated as 
this medium, confidence from any given model depends 
on its ability to describe multiple observables simultane- 
ously. The intermediate-mass range dilepton measure- 
ment ^9S\ at PHENIX is dominated by correlated decays 
of charm mesons which are not included in the calcu- 
lations of any published model. In the coming years, 
PHENIX will measure the yields and flow of bottom and 
charm mesons separately, which will further constrain the 
interpretation of the data. 



APPENDIX: DATA TABLES 

Data are presented in this Appendix for the minimum 
bias events (0-92%) and for each centrality class (0-10%, 
10-20%, 20-40%, 40-60%, and 60-92%) for Au-t-Au colli- 
sions at ^/suN = 200 GeV at midrapidity. Tables IXllXIII 
tabulate the differential invariant yield of heavy flavor 
electrons. Tables IXIIHIXVI give the nuclear modification 
factor -RauAu of heavy fiavor electrons. Table lXVlIl shows 
V2 for heavy flavor electrons. Table IXVIIII gives the dif- 
ferential invariant cross section of heavy flavor electrons. 
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TABLE X: DifFerential invariant yield of electrons {{N^+ +N^-)/2) from heavy-flavor decays for (upper, minimum bias) 0-92% 
and (lower) 0-10% centrality classes. The is in units of GeV/c. The yield and corresponding errors are in units of (GeV/c)~^ 



Centrality 



Pt 



Invariant yield 
7.30 X 10-^ 

3.74 X 10-^ 
2.09 X 10-^ 
1.22 X 10~^ 
7.46 X 10"^ 

5.04 X 10-=* 
3.27 X 10-=* 
1.84 X 10-=* 

8.21 X IQ-* 

3.90 X 10-* 
1.87 X 10-* 

9.49 X 10-^ 

5.05 X 10-^ 

2.91 X lO"'"^ 
1.71 X 10"'' 
9.98 X 10-"^ 

6.07 X 10-'^ 
3.87 X 10-** 
2.45 X 10-'^ 

1.61 X 10"'' 
1.25 X lO"** 
7.05 X lO-'' 
3.81 X lO-'' 
1.48 X lO-'' 
4.04 X 10-* 
1.18 X 10-"* 
5.84 X 10-" 
1.93 X 10-^ 
3.35 X 10-^ 

1.62 X 10-^ 

8.75 X 10-^ 
4.83 X 10-^ 
2.96 X 10-^ 
1.83 X 10-^ 
1.29 X 10-^ 
6.78 X 10-=* 
3.25 X 10-=* 

1.50 X 10"=* 
6.62 X 10-* 
3.32 X 10-* 
1.70 X 10-* 
9.52 X 10-^ 
5.45 X lO-'^ 
3.32 X 10-** 
2.03 X 10-^ 
1.11 X 10-^ 
7.73 X 10-^ 
4.96 X 10-*^ 
3.65 X 10-'^ 
2.56 X 10-^ 

1.08 X 10-** 

3.76 X lO-'' 

1.22 X 10"^ 
1.21 X 10"** 
2.38 X 10-* 
7.78 X 10-" 



Stat, error (+) 

3.01 X 10-* 
1.29 X 10-* 

6.55 X 10-* 
3.86 X 10"* 
2.32 X 10-* 
1.54 X 10-* 
1.06 X 10-* 
4.74 X 10-^ 
2.04 X lO-'' 
8.81 X 10-*^ 

2.56 X 10-^ 
1.35 X 10-^ 
7.71 X lO-'' 
4.83 X 10"^ 
3.12 X 10"^ 
2.12 X 10"^ 
1.50 X 10-'' 
1.11 X lO-'' 
8.27 X 10"* 
6.39 X 10"* 
5.23 X 10"* 
3.94 X 10-* 
1.71 X 10-* 
1.03 X 10-* 
4.45 X 10-" 
2.22 X 10-" 
1.59 X 10-" 

1.02 X 10-" 



Stat, error (— ) Sys. error (+) Sys. error ( — ) 



2.97 X 10-* 
1.27 X 10-* 
6.47 X 10-* 

3.82 X 10"* 
2.30 X 10-* 

1.52 X 10-* 
1.05 X 10-* 
4.60 X 10-^ 
1.99 X 10-^ 
8.62 X 10-*^ 
2.56 X 10-** 
1.35 X 10-*^ 
7.71 X lO-'' 

4.83 X 10-'' 
3.12 X 10-'' 
2.12 X 10"^ 
1.50 X 10-'' 
1.11 X 10-'' 
8.27 X 10-* 
6.39 X 10-* 
5.23 X 10"* 
3.90 X 10-* 
1.71 X 10-* 
1.02 X 10-* 
4.37 X 10-" 
2.15 X 10-" 

1.53 X 10-" 
6.87 X 10 



1.94 X 10-^ 
7.73 X 10-* 

3.66 X 10-* 

1.94 X 10-* 
1.04 X 10-* 
6.09 X 10-* 
3.70 X 10-* 

1.95 X 10-* 

9.04 X 10"'' 
4.35 X 10"" 
2.87 X 10-^ 
1.40 X 10-^ 
7.18 X 10-" 

3.90 X 10-" 
2.18 X 10-" 
1.26 X 10"" 
7.54 X 10-'' 

4.67 X 10-'' 
2.93 X 10"'" 

1.91 X 10"' 
1.38 X 10-'' 
8.40 X 10-* 
4.38 X 10-* 
1.76 X 10-* 
5.61 X 10-" 
1.61 X 10-" 
6.93 X 10-^° 
3.37 X 10-" 
6.70 X 10-^ 

2.73 X 10-^ 
1.30 X 10"^ 
6.87 X 10"* 

3.74 X 10-* 
2.17 X 10-* 

1.37 X 10-* 

7.05 X 10"* 
4.01 X 10"* 
1.91 X 10-* 
9.80 X 10-^ 

4.68 X 10-^ 
2.34 X 10-^ 
1.25 X 10-^ 
6.89 X 10"" 

4.06 X 10-" 
2.43 X 10-" 

1.38 X 10-" 
9.14 X 10-'' 

5.87 X 10"'" 
4.08 X 10"'' 
2.84 X 10-'' 
1.28 X 10-'' 

4.88 X 10-* 
1.72 X 10"* 
3.20 X 10"" 
2.66 X 10"" 
1.25 X 10-" 



1.91 X 10-^ 
7.63 X 10-* 
3.61 X 10-* 

1.92 X 10"* 
1.03 X 10-* 
6.03 X 10-* 
3.67 X 10-* 
1.90 X 10-* 
8.81 X lO"-"* 
4.26 X 10"^ 
2.87 X 10-^ 
1.40 X 10-^ 
7.18 X 10-" 

3.90 X 10"" 
2.18 X 10"" 
1.26 X 10"" 
7.54 X 10-'' 

4.67 X 10-'' 

2.93 X 10"'' 

1.91 X 10"^ 
1.38 X 10-'' 
8.40 X 10-* 
4.38 X 10-* 
1.76 X 10-* 
5.61 X 10-" 
1.61 X 10"" 
6.93 X 10"" 

3.37 X 10-" 
6.60 X 10-^ 

2.69 X 10-^ 
1.28 X 10-^ 

6.79 X 10-* 

3.70 X 10-* 
2.14 X 10-* 
1.36 X 10-* 

6.83 X 10-* 
3.90 X 10-* 

1.86 X 10-* 

9.80 X 10-^ 

4.68 X 10-^ 
2.34 X 10-^ 
1.25 X 10-^ 
6.89 X 10"" 
4.06 X 10-" 
2.43 X 10-" 

1.38 X 10-" 
9.14 X 10-'' 

5.87 X 10"'' 
4.08 X 10"'' 

2.84 X 10-'' 
1.28 X 10-'' 

4.88 X 10-* 
1.72 X 10"* 
3.20 X 10"" 
2.66 X 10"" 
1.25 X 10-" 



0-92% 



0.30 
0.50 
0.60 
0.70 
0.80 
0.80 
0.90 
1.10 
1.30 
1.50 
1.70 
1.90 
2.10 
2.30 
2.50 
2.70 
2.90 
3.10 
3.30 
3.50 
3.70 
3.90 
4.20 
4.80 
5.50 
6.50 
7.50 
8.50 



-10 



1.63 X 10-^ 
7.05 X 10-* 
3.59 X 10-* 
2.15 X 10-* 

1.31 X 10-* 
8.65 X 10-* 
6.13 X 10-* 
2.74 X 10"* 

7.08 X 10"^ 

3.09 X 10"^ 
1.05 X 10-** 
5.84 X 10-" 

3.51 X 10-" 
2.30 X 10-" 
1.55 X 10-" 
1.11 X 10-" 
8.02 X 10-'' 
5.84 X 10-'' 
4.53 X 10-'' 

3.52 X 10-'' 
2.82 X 10-'' 

2.32 X 10-'' 
9.32 X 10-* 
5.51 X 10-* 
2.50 X 10"* 
1.11 X 10"* 
1.15 X 10"* 
6.97 X 10-" 



1.61 X 10-^ 
6.95 X 10-* 
3.55 X 10-* 
2.12 X 10-* 
1.30 X 10-* 
8.57 X 10-* 
6.07 X 10-* 
2.66 X 10-* 
6.89 X 10-" 
3.02 X lO""^ 
1.05 X 10-^ 
5.84 X 10-" 
3.51 X 10-" 
2.30 X 10-" 
1.55 X 10-" 
1.11 X 10-" 
8.02 X 10-'' 
5.79 X 10-'' 
4.48 X 10-'' 
3.48 X 10-'' 
2.78 X 10-'' 
2.28 X 10-'' 
9.19 X 10-* 
5.39 X 10-* 
2.41 X 10-* 
8.26 X 10-" 
7.78 X 10-" 
3.70 X 10-" 



0-10% 



0.30 
0.50 
0.60 
0.70 
0.80 
0.80 
0.90 
1.10 
1.30 
1.50 
1.70 
1.90 
2.10 
2.30 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 
8.5 
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TABLE XI: Differential invariant yield of electrons {{N^+ +N^-)/2) from heavy-flavor decays for (upper) 10-20% and (lower) 
20-40% centrality classes. The px is in units of GeV/c. The yield and corresponding errors are in units of (GeV/c)~^ 



Centrality 



Pt 
0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.10 
2.3 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 
8.5 



Invariant yield 
1.90 X 10"^ 
9.11 X 10"^ 
5.29 X IQ-^ 

2.92 X IQ-^ 
1.88 X IQ-^ 
1.29 X 10"^ 

7.93 X 10"^ 
4.47 X IQ-^ 
2.20 X IQ-^ 
1.01 X IQ-^ 
4.69 X 10"" 

2.41 X 10"" 
1.28 X 10"" 
7.54 X 10"^ 
4.18 X 10-^ 
2.38 X 10"*^ 
1.40 X 10"^ 
9.37 X lO"'' 
5.59 X 10-*^ 
3.98 X 10-** 
3.14 X 10-*^ 
1.40 X 10"*^ 
9.77 X 10"^ 

3.42 X 10"'' 
9.71 X 10"* 
4.00 X 10-* 
2.45 X 10-* 
7.54 X 10-° 



Stat, error (+) Stat, error (— ) 



1.09 X 10-^ 
4.75 X 10-^ 
2.49 X 10-^ 
1.51 X 10-^ 
9.42 X 10-* 
6.36 X 10-* 
4.46 X 10-* 
2.01 X 10-* 

4.90 X 10-^ 
2.19 X 10-^ 
7.98 X 10-*^ 
4.54 X lO-'^ 
2.81 X 10"'' 

1.89 X 10-** 
1.27 X 10-** 
8.92 X 10-'' 
6.44 X lO-'' 

4.91 X 10"^ 
3.65 X 10-'' 

2.90 X 10-''' 
2.44 X 10-'' 
1.68 X 10-'' 

8.01 X 10-* 
4.70 X 10-* 

2.02 X 10-* 
1.13 X 10-* 
7.40 X 10-° 
3.30 X 10-° 



Sys. error (+) 
4.48 X 10-^ 

1.79 X 10-^ 
8.57 X 10-^ 
4.55 X 10-^ 
2.47 X 10-^ 

1.47 X 10"-'' 
8.77 X 10"* 

4.69 X 10-* 
2.50 X 10-* 
1.16 X 10-* 
6.73 X 10-^ 

3.27 X 10"^ 

1.67 X 10"-°^ 

9.24 X 10-" 

5.04 X 10-" 
2.86 X 10"" 
1.66 X 10"'' 

1.06 X 10"" 
6.41 X 10-'' 
4.39 X 10-'' 

3.25 X 10-'' 

1.68 X 10-'' 

1.05 X 10"^ 

3.92 X 10-* 

1.26 X 10-* 

4.53 X 10-° 
2.38 X 10-° 

1.07 X 10"° 

2.48 X 10"^ 
9.85 X 10-* 
4.66 X 10-^ 

2.50 X 10-^ 
1.33 X lO"-"* 

7.83 X 10"* 
4.66 X 10-* 

2.51 X 10-* 
1.26 X 10-* 
6.07 X 10-^ 
3.76 X 10-^ 
1.85 X 10-^ 
9.60 X 10-" 
5.15 X 10-" 

2.93 X 10-" 

1.70 X 10-" 
1.03 X 10"" 
6.43 X 10"^ 

4.09 X 10-'' 

2.54 X 10-'' 

1.80 X 10-'' 

1.10 X 10"'' 

5.84 X 10"** 
2.47 X 10"* 
7.65 X 10-° 

2.28 X 10-° 
5.98 X 10-1° 
3.98 X 10-1° 



Sys. error ( — " 

4.42 X 10-^ 
1.77 X 10-^ 
8.48 X 10-* 
4.50 X 10-* 

2.45 X 10-* 

1.46 X lO"-'' 
8.70 X 10-* 
4.55 X 10-* 

2.44 X 10-* 
1.13 X 10-* 

6.73 X 10-^ 

3.27 X 10-^ 

1.67 X 10"^ 

9.24 X 10-" 

5.04 X 10-" 
2.86 X 10"" 
1.66 X 10"" 

1.06 X 10"" 
6.41 X 10-'' 
4.39 X 10-'' 

3.25 X 10-'' 

1.68 X 10-'' 

1.05 X 10"'' 

3.92 X 10-* 

1.26 X 10-* 

4.53 X 10-° 
2.38 X 10-° 

1.07 X 10"° 

2.45 X 10"^ 

9.74 X 10-* 
4.61 X 10-* 
2.48 X 10-* 
1.31 X 10"* 
7.76 X 10"* 
4.63 X 10-* 
2.44 X 10-* 
1.23 X 10-* 
5.95 X 10-^ 
3.76 X 10-^ 
1.85 X 10-^ 
9.60 X 10-" 
5.15 X 10-" 

2.93 X 10-" 
1.70 X 10-" 
1.03 X 10"" 

6.43 X 10"'' 

4.09 X 10-'' 

2.54 X 10-'' 
1.80 X 10-'' 

1.10 X 10"'' 
5.84 X 10"* 

2.47 X 10"* 
7.65 X 10-° 

2.28 X 10-° 
5.98 X 10-1° 
3.98 X 10-1° 



10-20% 



1.11 X 10 
4.81 X 10-* 

2.52 X 10-* 

1.53 X 10-* 
9.51 X 10-* 
6.42 X 10-* 
4.50 X 10-* 
2.07 X 10-* 
5.03 X 10-^ 
2.24 X 10-^ 
7.98 X 10"" 

4.54 X 10"" 

2.81 X 10"" 
1.89 X 10-" 
1.27 X 10-" 
8.92 X 10"^ 
6.44 X 10"^ 
4.96 X 10-'' 
3.69 X 10-'' 
2.94 X 10-'' 
2.48 X 10-'' 
1.71 X 10-'' 
8.13 X 10-* 

4.82 X 10-* 
2.10 X 10-* 
1.20 X 10-* 
1.13 X 10-* 
6.78 X 10"° 



5.09 X 10-* 
2.28 X 10-* 
1.19 X 10-* 
7.47 X 10-* 
4.59 X 10-* 
3.13 X 10-* 
2.19 X 10-* 
9.98 X 10-** 
2.70 X 10-^ 
1.21 X 10-^ 
4.27 X 10-" 
2.42 X 10-" 

1.49 X 10-" 
9.78 X 10-'' 
6.70 X 10-'' 
4.73 X 10-'' 

3.50 X 10-'' 
2.64 X 10-'' 
2.00 X 10-'' 
1.52 X 10-'' 

1.25 X 10-'' 

9.26 X 10-* 
4.15 X 10-* 
2.52 X 10"* 
1.07 X 10-* 
5.45 X 10-° 
2.32 X 10-° 
1.36 X 10-° 



20-40% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.10 
2.3 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 
8.5 



7.35 X 10-^ 
4.19 X 10-^ 
2.25 X 10-^ 

1.55 X 10-^ 
8.93 X 10"* 
6.50 X 10"* 
3.79 X 10-* 
2.38 X 10-* 
1.15 X 10-* 
5.45 X 10-* 

2.56 X 10-* 
1.31 X 10-* 
7.12 X 10-" 
4.00 X 10-" 
2.35 X 10-" 
1.41 X 10-" 



8.86 : 



) X 10-" 

5.62 X 10-" 

3.63 X 10-" 
2.25 X 10-" 



1 



X 10-" 
9.65 X 10-^ 
5.37 X 10-^ 
2.27 X 10-^ 
6.08 X 10-* 
1.91 X 10-* 
3.95 X 10-° 
2.39 X 10-° 



5.15 X 10"* 
2.30 X 10-* 

1.20 X 10-* 
7.55 X 10-* 

4.63 X 10"* 

3.16 X 10"* 

2.21 X 10"* 
1.03 X 10-* 

2.77 X 10-" 
1.24 X 10-" 
4.27 X 10-" 
2.42 X 10-" 

1.49 X 10-" 

9.78 X 10-'' 
6.70 X 10-'' 
4.73 X 10-'' 

3.50 X 10-'' 

2.64 X 10-'' 
2.02 X 10-'' 
1.54 X 10-'' 
1.26 X 10-'' 
9.42 X 10-* 
4.21 X 10"* 
2.57 X 10-* 
1.11 X 10-* 
5.75 X 10-° 
4.02 X 10-° 
2.85 X 10-° 
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TABLE XII: Differential invariant yield of electrons {{N^+ +N^-)/2) from heavy-flavor decays for (upper) 40-60% and (lower) 
60-92% centrality classes. The px is in units of GeV/c. The yield and corresponding errors are in units of (GeV/c)~^ 



CciiUaliL^" 



I'T 

0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.10 
2.3 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 



Iii\'a.i iaiil 



Slal.. rrror (+) 
2.66 X 10"^ 
1.20 X 10"^ 

6.44 X 10~* 
4.00 X 10"* 
2.51 X IQ-* 
1.71 X IQ-* 
1.23 X 10"*' 

5.70 X 10"^ 
1.05 X 10"** 

5.00 X 10-^ 

2.01 X 10-^ 
1.22 X 10-^ 
7.91 X 10"'' 

5.45 X 10"'' 

3.85 X 10-'' 

2.71 X 10"'' 
1.98 X 10"'' 
1.60 X 10"'' 
1.18 X 10"'' 
9.34 X lO""* 

7.86 X 10"^ 
5.82 X 10"* 
2.62 X 10"** 
1.75 X lO""* 
7.04 X 10"" 
4.37 X lO"** 
2.74 X 10"^ 



til. 



I. error [ — ) 



2.63 X 10"^ 
1.18 X lO"'* 
6.38 X 10"" 
3.96 X 10"* 
2.49 X 10"* 

1.70 X 10"* 
1.22 X 10"* 
5.55 X lO"'"^ 
1.03 X 10"^ 

4.90 X 10"^ 
2.01 X 10"^ 
1.22 X 10"^ 

7.91 X 10"'' 
5.45 X 10"'' 
3.85 X 10"'' 

2.71 X 10"'' 
1.96 X 10"'' 
1.58 X 10"'' 
1.16 X 10"'' 
9.18 X 10"** 

7.70 X 10"* 
5.67 X 10"* 
2.57 X 10"* 

1.71 X 10"* 
6.71 X lO"** 
3.21 X 10"" 
1.13 X 10" 



S\',-.. tirr^ji' ( — 

9.42 X 10"^ 

3.69 X 10"^ 

1.72 X 10"^ 
9.02 X 10"* 

4.73 X 10"* 

2.74 X 10"* 
1.68 X 10"* 
8.65 X 10"^ 
4.23 X 10"^ 
2.10 X 10"^ 
1.33 X 10"^ 
6.67 X 10"^ 
3.54 X 10"" 
1.96 X 10"" 

1.13 X 10"" 
6.37 X 10"'' 

3.70 X 10"'' 
2.54 X 10"'' 

1.50 X 10"'' 
9.93 X 10"* 

7.43 X 10"* 

4.40 X 10"* 

2.41 X 10"* 
1.17 X 10"* 
3.21 X 10~" 
1.08 X 10"" 
2.04 X 10"" 
2.16 X 10"* 
8.06 X 10"* 
3.72 X 10"* 
1.93 X 10"* 
9.85 X 10"^ 

5.39 X 10"^ 
3.29 X 10"^ 
1.74 X 10"^ 
1.04 X 10"^ 
5.21 X 10"" 
2.27 X 10"" 

1.14 X 10"" 
5.54 X 10"'' 
3.16 X 10"'' 
1.89 X 10"'' 
1.08 X 10"'' 
6.46 X 10"* 

4.40 X 10"* 

2.58 X 10"* 
1.78 X 10"* 
1.33 X 10"* 
8.08 X 10"" 
3.82 X 10"" 
1.14 X 10"" 

3.59 X 10"" 

1.51 X 10"" 



9.31 X 10"-' 
3.65 X 10"-'' 
1.70 X 10"* 
8.94 X 10"* 

4.69 X 10"* 
2.72 X 10~* 
1.67 X 10"* 

8.42 X lO"'' 
4.13 X 10"^ 
2.06 X 10"^ 

1.33 X 10"^ 
6.67 X 10"" 
3.54 X 10"" 

1.96 X 10"" 

1.13 X 10"" 
6.37 X 10"'' 

3.70 X 10"'' 
2.54 X 10"'' 

1.50 X 10"'' 
9.93 X 10"* 

7.43 X 10"* 

4.40 X 10"* 

2.41 X 10"* 
1.17 X 10"* 
3.21 X 10"" 
1.08 X 10"" 
2.04 X 10"" 

2.14 X 10"* 

7.97 X 10"* 
3.69 X 10"* 
1.91 X 10"* 

9.77 X 10"^ 

5.34 X 10"^ 
3.27 X 10"^ 
1.69 X 10"^ 
1.01 X 10"^ 
5.10 X 10"" 
2.27 X 10"" 
1.14 X 10"" 
5.54 X 10"'' 
3.16 X 10"'' 
1.89 X 10"'' 
1.08 X 10"'' 
6.46 X 10"* 
4.40 X 10"* 

2.58 X 10"* 

1.78 X 10"* 
1.33 X 10"* 
8.08 X 10"" 
3.82 X 10"" 
1.14 X 10"" 

3.59 X 10"^" 

1.51 X 10"" 



40-60% 



2.13 X 10"^ 
1.36 X 10"^ 
8.48 X 10"* 
4.30 X 10"* 
2.78 X 10"* 
2.03 X 10"* 
1.41 X 10"* 
7.48 X 10"* 
4.05 X 10"* 
1.98 X 10"* 
8.78 X 10"^ 
4.64 X 10"^ 
2.59 X 10"^ 
1.50 X lO"'* 
9.10 X 10"" 
5.05 X 10"" 
2.90 X 10"" 

10"" 
10"" 
8.68 X 10"^ 

7.14 X 10"'' 
3.85 X 10"'' 
2.30 X 10"^ 
1.19 X 10"^ 
2.70 X 10"* 
1.03 X 10"* 
1.23 X 10"" 



2.27 
1.25 



8.47 X 10"* 
3.65 X 10"* 
1.99 X 10"* 
1.25 X 10"* 

7.81 X 10"^ 
5.24 X 10"" 
3.78 X 10"" 
1.74 X 10"" 
2.03 X 10"" 
1.06 X 10"" 
5.43 X 10"'' 
3.43 X 10"'' 
2.18 X 10"'' 
1.56 X 10"'' 
1.14 X 10"'' 
8.06 X 10"* 
5.97 X 10"* 
4.84 X 10"* 
3.52 X 10"* 
2.84 X 10"* 
2.43 X 10"* 
1.84 X 10"* 
7.63 X 10"" 
3.84 X 10"" 
1.37 X 10"" 

7.82 X 10"" 



60-92% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.10 
2.3 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 



5.28 X 10"* 
1.83 X 10"* 
1.18 X 10"* 

7.89 X 10"* 
4.27 X 10"* 
2.30 X 10"* 
1.99 X 10"* 
1.27 X 10"* 
6.15 X 10"" 
2.99 X 10"" 
1.39 X 10"" 
7.09 X 10"" 

3.08 X 10"" 
1.91 X 10"" 
1.25 X 10"" 
6.77 X 10"'' 

4.09 X 10"'' 
3.27 X 10"'' 
1.68 X 10"'' 
1.30 X 10"^ 
1.12 X 10"'' 
6.20 X 10"* 

2.90 X 10"* 
5.34 X 10"" 

1.29 X 10"" 
1.12 X 10"" 



8.56 X 10"* 
3.68 X 10"* 
2.01 X 10"* 
1.26 X 10"* 
7.87 X 10"" 
5.28 X 10"" 
3.81 X 10"" 
1.78 X 10"" 
2.08 X 10"" 
1.08 X 10"" 
5.43 X 10"'' 
3.43 X 10"'' 

2.18 X 10"'' 

1.57 X 10"^ 
1.15 X 10"^ 

8.19 X 10"* 
6.08 X 10"* 
4.95 X 10"* 
3.62 X 10"* 
2.94 X 10"* 
2.52 X 10"* 

1.92 X 10"* 

7.93 X 10"" 
4.10 X 10"" 
2.17 X 10"" 
1.48 X 10"" 
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TABLE Xlll: IiAn\n{l>'j ) fci' (iippci'r minimum bias) I) 1)2% and (Icnvcr) iO'/t cciiti'ality classes. 



Centrality 



Pt 



8.74 X 10-^ 
1.01 X 10+° 
1.01 X 10+° 

1.10 X 10+° 
9.34 X 10""- 

1.09 X 10+° 
9.49 X 10-^ 
9.44 X 10"^ 

1.06 X 10+° 
9.66 X 10"^ 

9.11 X 10-^ 
8.55 X 10-1 
8.00 X 10-1 
7.79 X 10-1 
7.47 X 10-1 
6.90 X 10-1 
6.47 X 10-1 
6.23 X 10-1 
5.84 X 10-1 

5.59 X 10-1 
6.18 X 10-1 
4.93 X 10-1 
4.74 X 10-1 
3.95 X 10-1 

3.07 X 10-1 

3.10 X 10-1 
4.68 X 10^1 
4.22 X 10-1 

1.08 X 10+° 
1.18 X 10+° 
1.14 X 10+° 

1.17 X 10+° 

1.00 X 10+° 

1.07 X 10+° 

1.01 X 10+° 
9.36 X 10-1 

1.13 X 10+° 
1.00 X 10+° 
8.71 X 10-1 

8.08 X 10-1 
7.29 X 10-1 

6.88 X 10-1 
6.42 X 10-1 

6.18 X 10-1 

5.83 X 10-1 

4.84 X 10-1 
4.97 X 10-1 
4.64 X 10-1 

4.89 X 10-1 
4.83 X 10-1 
3.62 X 10-1 
2.71 X 10-1 
2.49 X 10-1 
8.61 X 10-^ 

5.14 X 10-1 

4.60 X 10-1 



Stat, error (+) 



Stat, 



error (— ) 



Sys. error (+) 
3.19 X 10-1 
3.43 X 10-1 
2.65 X 10-1 
2.64 X 10-1 

1.64 X 10^1 

1.77 X 10-1 
1.30 X 10-1 

1.17 X 10-1 

1.35 X IQ-i 

1.16 X 10^1 

1.93 X 10-1 

1.78 X 10-1 
1.63 X 10-1 
1.53 X 10-1 

1.45 X 10-1 
1.30 X 10-1 
1.21 X 10-1 
1.15 X 10-1 
1.07 X 10-1 

1.02 X 10-1 
1.10 X 10-1 

9.03 X 10-^ 

8.55 X 10- = 
7.24 X 10- = 
6.03 X 10- = 

6.06 X 10- = 

8.56 X 10" = 
9.42 X 10- = 
4.37 X 10-1 
4.37 X 10-1 

3.18 X 10-1 
2.96 X 10-1 

1.80 X 10-1 
1.78 X 10-1 
1.42 X 10-1 

1.19 X 10-1 

1.36 X 10-1 
1.15 X 10-1 

1.81 X 10-1 

1.65 X 10-1 

1.46 X 10-1 
1.34 X 10-1 
1.24 X 10-1 
1.15 X IQ-i 

1.07 X 10-1 
9.01 X 10-= 
9.07 X 10- = 

8.47 X 10- = 
8.72 X 10" = 
8.60 X 10- = 
6.62 X 10- = 
5.15 X 10-= 

4.94 X 10- = 
2.58 X 10- = 

9.17 X 10- = 
9.77 X 10-= 



Sys. error ( — ) 
2.89 X 10-1 
3.15 X 10-1 
2.46 X 10"! 

2.46 X 10"! 
1.54 X 10"! 
1.68 X 10-1 
1.24 X 10-1 
1.12 X 10-1 
1.29 X 10-1 
1.11 X 10"! 

1.76 X 10-1 
1.62 X 10-1 

1.48 X 10-1 
1.38 X 10-1 

1.29 X 10"! 
1.15 X 10-1 
1.07 X 10-1 
1.01 X 10-1 
9.44 X 10- = 
8.99 X 10" = 

9.58 X 10" = 

7.96 X 10- = 

7.50 X 10- = 
6.38 X 10- = 
5.41 X 10- = 

5.43 X 10" = 
7.53 X 10" = 

8.68 X 10- = 
4.03 X 10-1 
4.07 X 10-1 

2.98 X 10-1 

2.77 X 10"! 
1.70 X 10-1 

1.69 X 10-1 
1.35 X 10-1 
1.14 X 10-1 

1.30 X 10"! 

1.10 X 10-1 
1.65 X 10-1 

1.49 X 10-1 
1.32 X 10-1 
1.20 X 10-1 

1.11 X 10-1 

1.01 X 10-1 

9.44 X 10- = 

7.97 X 10- = 

7.99 X 10- = 

7.45 X 10- = 
7.62 X 10" = 

7.51 X 10- = 
5.83 X 10- = 

4.59 X 10- = 
4.44 X 10- = 

2.47 X 10" = 

8.02 X 10" = 
8.92 X 10- = 



0-92% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.1 
2.3 
2.50 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 
8.5 



4.27 X 10-1 
3.96 X 10-1 

2.96 X 10-1 
3.12 X 10-1 
2.01 X 10-1 
2.42 X 10-1 
1.87 X 10-1 
1.29 X 10-1 
2.15 X 10-1 
2.27 X 10"! 
1.25 X 10- = 

1.21 X 10- = 

1.22 X 10- = 
1.29 X 10- = 
1.36 X 10- = 
1.46 X 10- = 
1.60 X 10- = 
1.79 X 10- = 

1.97 X 10- = 
2.21 X 10- = 
2.59 X 10- = 
2.76 X 10" = 
2.12 X 10- = 
2.75 X 10- = 
3.38 X 10- = 
5.85 X 10- = 
1.27 X 10"! 

2.23 X 10-1 



2.18 
2.24 
1.89 
2.01 
1.42 
1.70 
1.37 
1.03 
1.54 
1.55 
1.25 
1.21 
1.22 
1.29 
1.36 
1.46 
1.60 
1.79 
1.97 
2.21 
2.59 
2.73 
2.12 
2.71 
3.32 
5.67 
1.22 
1.50 



X 10-1 
X 10-1 

X 10"! 
X 10"! 
X 10"! 

X 10-1 
X 10-1 

X 10"! 
X 10"! 
X 10"! 

X 10-= 
X 10-= 
X 10-= 
X 10-= 

X 10" = 
X 10- = 
X 10-= 
X 10-= 
X 10-= 
X 10" = 
X 10" = 
X 10" = 
X 10-= 
X 10-= 
X 10-= 
X 10-= 
X 10"! 
X 10-1 



0-10% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.1 
2.30 
2.5 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 
8.5 



5.30 X 10-1 

4.65 X 10-1 
3.36 X 10-1 
3.36 X 10-1 
2.17 X 10-1 
2.40 X 10-1 
2.03 X 10-1 
1.32 X 10-1 
2.29 X 10-1 
2.35 X 10-1 

1.38 X 10-= 
1.42 X 10-= 
1.50 X 10- = 

1.66 X 10- = 
1.83 X 10- = 
2.07 X 10- = 

2.31 X 10- = 
2.54 X 10- = 
2.91 X 10- = 
3.29 X 10- = 
3.77 X 10" = 

4.39 X 10- = 
3.13 X 10- = 
3.97 X 10- = 
5.13 X 10- = 
7.88 X 10- = 
2.49 X 10-1 
4.12 X 10-1 



2.72 
2.64 
2.15 
2.18 
1.55 
1.71 
1.50 
1.06 
1.64 
1.61 
1.38 
1.42 
1.50 
1.66 
1.83 
2.07 
2.31 
2.51 
2.88 
3.25 
3.72 
4.31 
3.08 
3.88 
4.95 
5.88 
1.68 
2.18 



X 10-1 
X 10-1 
X 10-1 
X 10-1 
X 10-1 
X 10-1 
X 10-1 
X 10-1 

X 10"! 
X 10"! 

X 10-= 
X 10-= 
X 10-= 

X 10-= 
X 10" = 
X 10" = 
X 10-= 
X 10-= 
X 10-= 
X 10-= 
X 10" = 
X 10-= 
X 10-= 
X 10-= 
X 10-= 
X 10" = 
X 10-1 
X 10-1 
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TABLE XIV: ii'AuAu(i)j ) for (upper) 10 2()'/{ and (low(-r) 20 lO'X coutrality classct 



Centrality 



Pt 




Stat. 


error (+) 


Stat, error (- 


0.3 


9.72 X 


10-1 


4.76 


X 


10-1 


2.46 X 10-1 


0.5 


1.05 X 


10+° 


4.13 


X 


10-1 


2.36 X 10-1 


0.6 


1.09 X 


10+° 


3.22 


X 


10-1 


2.08 X IQ-i 


0.7 


1.12 X 


10+° 


3.22 


X 


10-1 


2.11 X 10-1 


0.8 


1.01 X 


10+° 


2.20 


X 


10-1 


1.58 X 10-1 


0.8 


1.19 X 


10+° 


2.67 


X 


10-1 


1.91 X 10-1 


0.9 


9.82 X 


10-1 


1.99 


X 


10-1 


1.48 X 10-1 


1.1 


9.76 X 


10-1 


1.39 


X 


10-1 


1.12 X 10-1 


1.3 


1.22 X 


10+° 


2.45 


X 


10-1 


1.76 X IQ-i 


1.5 


1.07 X 


10+" 


2.51 


X 


10-1 


1.72 X IQ-i 


1.7 


9.74 X 


10-1 


1.66 


X 


10-^ 


1.66 X 10-^ 


1.9 


9.25 X 


10-1 


1.75 


X 


10-^ 


1.75 X 10-^ 


2.1 


8.68 X 


10-1 


1.90 


X 


10-^ 


1.90 X 10-^ 


2.3 


8.61 X 


10-1 


2.15 


X 


10-^ 


2.15 X 10-^ 


2.50 


7.78 X 


10-1 


2.37 


X 


10-^ 


2.37 X 10-^ 


2.7 


7.00 X 


10-1 


2.63 


X 


10-^ 


2.63 X 10-^ 


2.9 


6.39 X 


10-1 


2.93 


X 


10-^ 


2.93 X 10-^ 


3.1 


6.43 X 


10-1 


3.40 


X 


10-^ 


3.37 X 10-^ 


3.3 


5.67 X 


10-1 


3.74 


X 


10-^ 


3.70 X 10-^ 


3.5 


5.88 X 


10-1 


4.34 


X 


10-^ 


4.29 X 10-^ 


3.7 


6.65 X 


10-1 


5.23 


X 


10-^ 


5.16 X IQ-^ 


3.9 


4.19 X 


10-1 


5.11 


X 


10-^ 


5.01 X IQ-^ 


4.2 


5.18 X 


10-1 


4.31 


X 


10-^ 


4.24 X 10-^ 


4.8 


3.89 X 


10-1 


5.48 


X 


10-^ 


5.35 X 10-^ 


5.5 


3.14 X 


10-1 


6.79 


X 


10-^ 


6.54 X 10-^ 


6.5 


4.50 X 


10-1 


1.35 


X 


10-1 


1.27 X 10-1 


7.5 


8.39 X 


10-1 


3.88 


X 


IQ-i 


2.53 X 10-1 


8.5 


7.03 X 


10-1 


6.32 


X 


10-1 


3.08 X 10-1 


0.3 


7.65 X 


10-1 


3.76 


X 


10-1 


1.96 X 10-1 


0.5 


9.79 X 


10-1 


3.87 


X 


10-1 


2.21 X 10-1 


0.6 


9.47 X 


10-1 


2.80 


X 


10-1 


1.81 X 10-1 


0.7 


1.21 X 


10+° 


3.48 


X 


10-1 


2.27 X 10-1 


0.8 


9.71 X 


10-1 


2.12 


X 


10-1 


1.53 X 10-1 


0.8 


1.23 X 


10+° 


2.75 


X 


10-1 


1.96 X 10-1 


0.9 


9.56 X 


10-1 


1.94 


X 


10-1 


1.45 X 10-1 


1.1 


1.06 X 


10+° 


1.49 


X 


10-1 


1.21 X 10-1 


1.3 


1.29 X 


10+° 


2.61 


X 


10-1 


1.87 X IQ-i 


1.5 


1.17 X 


10+° 


2.75 


X 


10-1 


1.88 X 10-1 


1.7 


1.08 X 


10+° 


1.81 


X 


10-^ 


1.81 X 10-^ 


1.9 


1.03 X 


10+° 


1.89 


X 


10-^ 


1.89 X 10-^ 


2.1 


9.81 X 


10-1 


2.05 


X 


10-^ 


2.05 X 10-^ 


2.3 


9.30 X 


10-1 


2.27 


X 


10-^ 


2.27 X 10-^ 


2.5 


8.91 X 


10-1 


2.54 


X 


10-^ 


2.54 X 10-^ 


2.7 


8.46 X 


10-1 


2.84 


X 


10-^ 


2.84 X 10-^ 


2.9 


8.20 X 


10-1 


3.24 


X 


10-^ 


3.24 X 10-^ 


3.1 


7.86 X 


10-1 


3.68 


X 


10-^ 


3.68 X 10-^ 


3.3 


( .OU X 


in-i 

J-U 


4.17 


X 


1 


4.13 X 10-^ 


3.5 


6.77 X 


10-1 


4.63 


X 


10-^ 


4.57 X 10-^ 


3.7 


7.26 X 


10-1 


5.44 


X 


10-^ 


5.37 X 10-^ 


3.9 


5.87 X 


10-1 


5.73 


X 


10-^ 


5.63 X 10-^ 


4.2 


5.80 X 


10-1 


4.54 


X 


10-^ 


4.48 X 10-^ 


4.8 


5.26 X 


10-1 


5.96 


X 


10-^ 


5.84 X 10-^ 


5.5 


4.01 X 


10-1 


7.31 


X 


10-^ 


7.08 X 10-^ 


6.5 


4.36 X 


10-1 


1.32 


X 


10-1 


1.25 X lO"! 


7.5 


2.75 X 


10-1 


2.80 


X 


10-1 


1.62 X 10-1 


8.5 


4.54 X 


10-1 


5.42 


X 


10-1 


2.58 X 10-1 



Sys. error (+) 
3.66 X 10-1 

3.62 X 10-1 
2.92 X 10-1 
2.71 X 10-1 
1.77 X 10-1 
1.95 X 10-1 
1.34 X 10-1 
1.21 X 10-1 

1.48 X 10-1 
1.21 X IQ-i 

2.00 X 10-1 
1.85 X 10-1 
1.70 X 10-1 

1.63 X 10-1 
1.47 X IQ-i 
1.29 X 10-1 
1.17 X 10-1 
1.15 X 10-1 
1.02 X 10-1 

1.04 X IQ-i 
1.15 X 10-1 
7.68 X 10-^ 
9.10 X 10-^ 
7.02 X 10-^ 
5.99 X 10-^ 
8.07 X 10-^ 

1.42 X IQ-i 

1.40 X 10-1 
2.82 X 10-1 
3.21 X 10-1 

2.43 X 10-1 
2.89 X IQ-i 

1.70 X 10-1 
1.97 X 10-1 

1.32 X 10-1 
1.29 X 10-1 

1.51 X 10-1 

1.33 X 10-1 
2.25 X 10-1 
2.09 X 10-1 
1.95 X 10-1 
1.80 X 10-1 

1.71 X 10-1 
1.55 X 10-1 

1.49 X 10-1 

1.41 X 10-1 

1.34 X 10-1 
1.21 X 10-1 

1.27 X lO"! 

1.05 X 10-1 
1.02 X 10-1 

9.28 X 10-^ 

7.52 X 10-^ 

8.01 X 10-^ 
5.65 X 10-^ 
9.85 X 10-^ 



Sys. error ( — ) 
3.34 X 10-1 
3.33 X 10-1 
2.72 X 10"! 
2.53 X IQ-i 
1.67 X 10"! 

1.85 X 10-1 
1.28 X 10-1 
1.16 X 10-1 
1.41 X IQ-i 

1.16 X IQ-i 
1.81 X 10-1 
1.67 X 10-1 
1.52 X 10-1 
1.45 X 10-1 

1.30 X IQ-i 

1.14 X 10-1 

1.03 X 10-1 
1.01 X 10-1 

8.97 X 10-^ 
9.11 X IQ-^ 

9.98 X 10-^ 
6.77 X 10-^ 

7.91 X 10-^ 

6.15 X 10-^ 
5.33 X 10-^ 
7.07 X IQ-^ 
1.22 X IQ-i 
1.26 X 10-1 
2.56 X 10-1 
2.93 X 10-1 
2.24 X IQ-i 
2.69 X 10-1 

1.60 X IQ-i 

1.86 X 10-1 

1.26 X 10-1 
1.24 X 10-1 
1.43 X 10-1 

1.27 X 10-1 

2.04 X 10-1 
1.90 X 10-1 
1.75 X 10-1 

1.61 X 10-1 
1.52 X 10-1 
1.37 X 10-1 

1.31 X 10-1 
1.24 X IQ-i 

1.17 X 10-1 
1.06 X 10-1 
1.11 X IQ-i 
9.24 X 10-^ 

8.92 X 10-^ 
8.09 X 10-^ 
6.67 X 10-^ 
7.06 X 10-^ 
5.13 X 10"^ 
9.03 X 10-^ 



10-20% 



20-40% 
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TABLE X\': ./?AuAu(i/j ) for (uiipcr) 10 (iOSt and flowc-i') (iO 927, ct-utralitv classes. 



Centrality 



Pt 



AAuAu 

7.27 X 10"^ 
1.04 X 10+° 
1.17 X 10+° 

1.10 X 10+° 
9.89 X 10-^ 
1.25 X 10+° 
1.17 X 10+° 
1.09 X 10+° 
1.49 X 10+° 
1.40 X 10+° 
1.22 X 10+° 
1.19 X 10+" 
1.17 X 10+° 
1.14 X 10+° 

1.13 X 10+° 
9.93 X 10"^ 

8.81 X 10""- 
1.04 X 10+° 

8.49 X 10-^ 

8.54 X 10"^ 
1.01 X 10+° 
7.66 X 10-^ 
8.12 X 10"^ 
9.03 X 10"^ 

5.82 X 10"^ 

7.75 X 10"^ 
2.79 X 10"^ 

1.11 X 10+° 
8.62 X 10-1 
1.01 X 10+° 
1.25 X 10+° 
9.38 X IQ-i 
8.77 X IQ-i 
1.01 X 10+" 

1.14 X 10+° 
1.40 X 10+° 
1.30 X 10+° 

1.19 X 10+° 

1.12 X 10+° 
8.56 X IQ-i 
8.95 X 10-1 

9.55 X 10-1 
8.21 X IQ-i 
7.66 X IQ-i 
9.24 X IQ-i 
7.01 X IQ-i 
7.92 X 10-1 

9.76 X 10-1 
7.62 X IQ-i 
6.33 X IQ-i 

2.50 X 10^1 
1.72 X 10-1 

5.20 X 10-1 



Stat, error (+) 



Stat 



error ( — ) 



Sys. error (+) 



3.23 X 10-1 
3.39 X 10-1 
2.99 X 10-1 
2.62 X 10-1 
1.80 X IQ-i 
2.01 X IQ-i 

1.59 X 10-1 
1.38 X 10-1 
1.83 X 10-1 

1.53 X IQ-i 
2.57 X IQ-i 
2.45 X 10^1 
2.33 X IQ-i 
2.22 X IQ-i 
2.16 X 10-1 
1.86 X 10-1 

1.66 X 10^1 

1.85 X IQ-i 
1.56 X 10-1 

1.54 X 10-1 
1.75 X 10-1 
1.38 X IQ-i 
1.41 X IQ-i 
1.54 X 10-1 
1.07 X 10-1 
1.35 X 10-1 
6.06 X IQ-^ 
4.14 X 10-1 
3.25 X 10-1 

2.78 X 10-1 
2.95 X 10-1 

1.86 X IQ-i 
1.74 X IQ-i 
1.49 X IQ-i 
1.44 X 10-1 
1.85 X 10-1 

1.67 X 10-1 

2.60 X 10-1 

2.44 X IQ-i 
1.98 X IQ-i 
1.97 X 10-1 

2.00 X 10-1 
1.73 X 10-1 

1.61 X IQ-i 

1.79 X IQ-i 

1.45 X IQ-i 
1.54 X 10-1 
1.78 X 10-1 
1.45 X IQ-i 
1.21 X IQ-i 
6.38 X 10"^ 
5.35 X 10-^ 

1.01 X 10-1 



Sys. error ( — 

3.03 X 10-1 

3.09 X 10-1 
2.76 X 10-1 
2.44 X 10-1 
1.71 X 10^1 
1.90 X 10-1 
1.51 X 10-1 
1.32 X 10-1 

1.75 X 10-1 

1.46 X 10-1 

2.34 X 10^1 

2.22 X 10^1 

2.10 X 10-1 
1.99 X 10-1 
1.92 X 10^1 
1.65 X 10^1 
1.48 X 10^1 
1.62 X 10-1 

1.37 X 10-1 

1.35 X 10-1 
1.51 X 10-1 
1.21 X 10-1 

1.23 X 10-1 
1.32 X 10-1 
9.39 X 10-^ 
1.17 X 10-1 

5.56 X 10''' 

3.76 X 10-1 

3.04 X 10-1 

2.59 X 10-1 
2.75 X 10-1 

1.77 X 10-1 
1.68 X 10-1 
1.43 X 10-1 

1.38 X 10-1 

1.78 X 10-1 
1.61 X 10-1 

2.39 X 10-1 

2.24 X 10^1 
1.84 X 10-1 

1.81 X 10-1 

1.82 X 10-1 
1.58 X 10-1 

1.47 X 10^1 

1.60 X 10-1 
1.32 X 10-1 
1.38 X 10-1 

1.57 X 10-1 
1.29 X 10^1 
1.08 X 10^1 
6.00 X 10-^ 
5.14 X 10-^ 
8.99 X 10-^ 



40-60% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.1 
2.3 
2.50 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 
7.5 



3.65 X 10-1 

4.18 X 10-1 

3.51 X 10-1 
3.28 X 10-1 
2.28 X 10-1 
2.94 X 10-1 
2.49 X 10-1 
1.68 X 10-1 
3.02 X 10-1 
3.28 X 10-1 
2.79 X 10-^ 
3.11 X 10"^ 
3.57 X 10-^ 

4.15 X 10-^ 
4.78 X 10-^ 
5.34 X 10-^ 
6.01 X 10"^ 

7.32 X 10-^ 
8.00 X 10-^ 

9.19 X 10-^ 
1.11 X 10-1 

1.16 X 10-1 

9.27 X 10-^ 

1.33 X 10-1 

1.52 X 10-1 

3.28 X 10-1 
6.24 X 10-1 



2.01 
2.46 
2.32 
2.23 
1.72 
2.18 
1.92 
1.41 
2.17 
2.25 
2.79 
3.11 
3.57 
4.15 
4.78 
5.34 
5.95 
7.24 
7.88 
9.03 
1.09 
1.13 
9.09 
1.30 
1.45 
2.40 
2.57 



X 10-1 
X 10-1 
X 10-1 
X 10-1 

X 10^1 

X 10"! 

X 10-1 
X 10-1 
X 10-1 

X IQ-i 
X 10"^ 
X 10"^ 
X 10-^ 
X 10-^ 
X 10-^ 
X 10-^ 
X 10"^ 
X 10-^ 
X 10-^ 
X 10-^ 
X 10-1 
X 10^1 
X 10"^ 
X IQ-i 
X IQ-i 
X 10-1 
X IQ-i 



60-92% 



0.3 
0.5 
0.6 
0.7 
0.8 
0.8 
0.9 
1.1 
1.3 
1.5 
1.7 
1.9 
2.1 
2.3 
2.50 
2.7 
2.9 
3.1 
3.3 
3.5 
3.7 
3.9 
4.2 
4.8 
5.5 
6.5 



5.69 X 10-1 
3.80 X 10-1 

3.39 X 10-1 

4.05 X 10-1 
2.64 X 10-1 
2.78 X 10-1 

2.77 X 10^1 
2.22 X 10-1 
2.84 X 10-1 
3.07 X 10-1 
4.64 X 10-^ 
5.43 X 10-^ 

6.06 X 10-^ 
7.38 X 10-^ 
8.84 X 10-^ 
9.93 X 10-^ 
1.14 X IQ-i 

1.40 X 10-1 
1.51 X 10-1 

1.78 X 10-1 
2.19 X 10-1 
2.36 X 10-1 
1.73 X 10-1 
1.92 X 10^1 
2.89 X 10-1 
6.86 X 10-1 



3.26 
2.56 
2.51 
3.00 
2.21 
2.40 
2.39 
1.97 
2.05 
2.12 
4.64 
5.43 
6.06 
7.31 
8.73 
9.78 
1.12 
1.37 
1.47 
1.73 
2.11 
2.26 
1.66 
1.80 
1.83 
3.62 



X 10-1 
X 10-1 
X IQ-i 
X 10-1 
X IQ-i 
X 10^1 
X 10^1 
X 10-1 
X 10-1 
X 10-1 
X 10-^ 
X 10"^ 
X 10-^ 
X 10-^ 
X 10-^ 
X 10-^ 
X 10^1 
X IQ-i 
X 10-1 
X 10-1 
X 10-1 
X 10^1 
X 10^1 
X 10^1 
X 10-1 
X 10-1 
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TABLE XVI: py-integrated nuclear modification factors -RauAu- 



Pt 




-RauAu 


Stat, error (+) 


Stat error ( — ) 


Sys. error (it) 


Common fractional error 




(GeV/c) 












from p+p cross section. 






14.5 


1.04 X 10+° 


7.64 X 10-^ 


7.55 X 10-^ 


3.88 X 10-1 








59.95 


9.86 X 10"^ 


3.90 X 10-^ 


3.86 X 10-^ 


2.53 X 10-1 






> 0.30 


140.4 


9.43 X 10"^ 


2.29 X 10-^ 


2.26 X 10-^ 


2.19 X 10-1 


+2.325 X 10-1 - 1.587 X 10" 






234.6 


1.03 X 10+° 


2.38 X 10-^ 


2.35 X 10-^ 


2.43 X 10-1 








325.2 


1.09 X 10+" 


2.20 X 10-^ 


2.17 X 10-^ 


2.71 X 10-1 








14.5 


1.08 X 10+" 


7.16 X 10-^ 


7.09 X 10-^ 


3.49 X 10-1 








59.95 


1.14 X 10+° 


3.71 X 10-^ 


3.67 X 10-^ 


2.08 X 10-1 






> 0.60 


140.4 


1.10 X 10+° 


2.11 X 10-^ 


2.09 X 10-^ 


1.71 X 10-1 


+ 1.309 X 10-1 - 1.037 X 10" 






234.6 


1.05 X 10+° 


2.11 X 10-^ 


2.08 X 10-^ 


1.59 X 10-1 








325.2 


1.03 X 10+" 


1.84 X 10-^ 


1.81 X 10-^ 


1.59 X IQ-i 








14.5 


1.07 X 10+" 


7.27 X 10-^ 


7.17 X 10-^ 


3.33 X 10-1 








59.95 


1.20 X 10+" 


3.82 X 10-^ 


3.77 X 10-^ 


1.97 X 10-1 






> 0.8 


140.4 


1.09 X 10+° 


2.14 X 10-^ 


2.11 X 10-^ 


1.49 X 10-1 


+1.096 X 10-1 - 8.989 X 10" 


-2 




234.6 


1.04 X 10+° 


2.13 X 10-^ 


2.09 X 10-^ 


1.35 X 10-1 








325.2 


9.77 X 10"^ 


1.81 X 10-^ 


1.78 X 10-^ 


1.29 X 10-1 








14.5 


7.78 X 10-^ 


2.96 X 10-^ 


2.92 X 10-^ 


2.57 X 10-1 








59.95 


9.91 X 10-^ 


1.66 X 10-^ 


1.66 X 10-^ 


1.77 X 10-1 






> 2.0 


140.4 


8.13 X 10-^ 


9.13 X 10-^ 


9.12 X 10-^ 


1.25 X 10-1 


+ 1.096 X 10-1 - 8.989 X 10" 


-2 




234.6 


7.10 X 10-^ 


8.52 X 10-^ 


8.50 X 10-^ 


1.03 X 10-1 








325.2 


5.89 X 10"^ 


6.64 X lO"-"* 


6.63 X 10"^ 


8.88 X 10-^ 








14.5 


7.04 X 10^^ 


6.39 X 10"^ 


5.98 X 10"^ 


2.24 X 10-1 








59.95 


8.55 X 10"^ 


3.35 X 10-^ 


3.25 X 10-^ 


1.40 X 10-1 






> 3.0 


140.4 


6.52 X 10"^ 


1.67 X 10-^ 


1.64 X 10-^ 


9.09 X 10-^ 


+9.049 X 10"^ - 7.662 x 10" 


-2 




234.6 


5.40 X 10-1 


1.56 X 10-^ 


1.52 X 10-^ 


7.12 X 10-^ 








325.2 


4.24 X 10^1 


1.17 X 10-^ 


1.14 X IQ-^ 


5.90 X 10-^ 








14.5 


3.99 X 10^1 


1.22 X 10^1 


9.72 X 10"^ 


1.31 X 10-1 








59.95 


7.59 X 10"^ 


6.73 X 10"^ 


6.17 X 10"^ 


1.22 X 10^1 






> 4.0 


140.4 


5.11 X 10-1 


3.16 X 10-^ 


3.00 X 10-^ 


7.14 X 10-^ 


+0.066 X 10-1 - 8.787 X 10" 


-2 




234.6 


4.51 X 10-1 


3.09 X 10-^ 


2.88 X 10-^ 


5.96 X 10-^ 








325.2 


3.03 X 10-1 


2.19 X 10-^ 


2.05 X 10-^ 


4.42 X 10-^ 
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TABLE XVII: Heavy-flavor V2 from Au+Au collisions, for 
the centralities indicated. The pr is in units of GeV/c. 



Centrality 


Pt 


V2 


stat error 


syst error 




0.546 


0.0401 


0.00636 


0.0189 




0.646 


0.0408 


0.00712 


0.0188 




0.746 


0.0385 


0.00742 


0.0168 




0.847 


0.0626 


0.00767 


0.0145 




0.947 


0.0625 


0.00861 


0.0137 




1.09 


0.0563 


0.00699 


0.0118 


0-92% 


1.29 


0.0698 


0.00923 


0.0115 




1.52 


0.092 


0.0103 


0.0112 




1.83 


0.087 


0.0152 


0.0107 




2.20 


0.0692 


0.0181 


0.00931 




2.70 


0.0706 


0.0308 


0.00842 




.3.24 


0.0308 


0.0465 


0.00789 




4.05 


0.00986 


0.0645 


0.00957 



TABLE XVIII: Differential invariant cross section of electrons 
((iVj,+ + iVg- ) /2) from heavy-flavor decays for 200 GeV p + p 
collisions at midrapidity. The px is in units of GeV/c. The 
cross section and corresponding errors are in units of mb 



Pt 


Invariant yield 


Stat, error 


Sys. error 


0.350 


1.36x10-^ 


4.45x10-^ 


5.95x10-^ 


0.400 


6.05x10-^ 


1.70x10-^ 


2.39x10-^ 


0.550 


3.36x10-^ 


7.56x10-* 


1.04x10-^ 


0.650 


1.81x10-^ 


3.98x10-* 


5.12x10-* 


0.750 


1.30x10"^ 


2.28x10-* 


2.59x10"* 


0.850 


7.50x10"" 


1.35x10-* 


1.35x10"* 


0.950 


5.61x10"* 


9.14x10"-'' 


7.85x10"^ 


1.10 


3.18x10-* 


3.77x10-*^ 


3.64x10-'' 


1.30 


1.26x10-* 


2.10x10-^ 


1.50x10-^ 


1.50 


6.58x10-^ 


1.25x10-^ 


6.56x10-" 


1.70 


3.97x10-^ 


2.07x10-'' 


3.32x10-" 


1.90 


1.99x10"-'^ 


1.18x10"'^ 


1.65x10"" 


2.10 


1.14x10-*^ 


7.44x10-'' 


8.96x10-'' 


2.30 


6.83x10-'' 


4.96x10-'' 


5.09x10-'' 


2.50 


3.98x10-** 


3.63x10-'' 


3.05x10-'' 


2.70 


2.44x10-** 


6.49x10-® 


2.80x10-'' 


2.90 


1.63x10-® 


4.77x10-® 


1.77x10"'' 


3.10 


1.05x10-'' 


3.53x10-® 


1.10x10"'' 


3.30 


7.21x10-'' 


2.73x10-® 


7.31x10-® 


3.50 


5.04x10-'' 


2.15x10-® 


4.89x10-® 


3.70 


3.45x10-'' 


1.70x10-® 


3.30x10-® 


3.90 


2.37x10-'' 


1.36x10-® 


2.25x10"® 


4.25 


1.33x10-'' 


5.95x10-^ 


1.22x10"® 


4.75 


6.13x10"** 


3.73x10"^ 


5.45x10"^ 


5.50 


2.06x10-** 


1.79x10-^ 


1.85x10-^ 


6.50 


6.62x10-^ 


9.26x10-1° 


5.72x10-1° 


7.50 


1.70x10-^ 


5.73x10"^" 


1.60x10"!" 


8.50 


9.72x10-^° 


5.16x10"" 


8.07xl0"ii 



0-10% 



0.40 
0.60 
0.80 
1.05 
1.40 
1.80 
2.50 
4.00 



5.13 
3.52 
4.13 

4.20 
4.08 
7.05 

9.81 
1.09 



xlQ-^ 
xlO-^ 
xlQ-^ 

xlO"^ 
xlO-^ 
xlO-^* 

xlO"-'' 

X 10 ^ 



1.16 
8.82 
9.13 

9.65 
1.02 
1.67 

2.22 
5.<-<3 



xlO"^ 
xlQ-^ 
xlQ-^ 
xlO"^ 
xlO-^ 
xlO-= 
xlO"^ 
X 1(1 - 



9.44 
6.77 
8.52 
5.45 
4.96 
7.54 
6.51 
l.:!l 



xlO"^ 
xlO-'* 
xlO-® 
xlO"^ 
xlO"^ 
xlO-® 
xlO"^ 
xll) - 
xlO"^ 
xlO"^ 
xlO"^ 
xlO-'* 
xlO-^ 
xlO"^ 
xlO"^ 
xlO-^ 
xlO-^ 
xlQ-^ 
xlQ-^ 
xlO-^ 
xlO"^ 
xlO"^ 
xlO"^ 
xlO-" 
xlO-^ 
xlO"^ 
xlO"^ 
xlO-^ 
xlO-^ 
xlO-^ 
xlO"^ 
xlO"^ 



10-20% 



0.40 
0.60 
0.80 
1.05 
1.40 
1.80 
2.50 
4.00 



2.90 
5.01 
4.78 
6.28 
5.89 
8.25 
9.60 
4.04 



xlO"^ 
xlO"^ 
xlO"2 
xlO-^ 
xlQ-^ 
xlO"^ 
xlO"^ 
xlO"^ 



9.59 
6.52 
6.59 
7.04 
7.65 
1.25 
1.63 
4.06 



xlO"-' 
xlO"-^ 
xlO"-'' 
xlQ-^ 
xlQ-^* 
xlO"^ 
xlO"^ 
xlO"^ 



2.03 
1.13 
1.07 
8.54 
8.19 
7.07 
6.34 
5.44 



20-40% 



0.40 
0.60 
0.80 
1.05 
1.40 
1.80 
2.50 
4.00 



6.55 
7.34 
1.12 
9.47 
1.30 
1.27 
8.48 
8.15 



xlQ-^ 
xlQ-^ 
xlQ-i 
xlO-^ 
xlO-i 
xlO-i 
xlQ-^ 
xlQ-^ 



8.01 
5.44 
5.73 
6.30 
7.08 
1.14 
1.47 
3.62 



xlQ-^ 
xlQ-^* 
xlQ-^* 
xlO-^* 

xlO"^ 
xlO"^ 
xlO"^ 
xlQ-'' 



1.82 
1.36 
1.05 
1.01 

8.69 
8.36 
8.17 
6.29 



40-60% 



0.40 
0.60 
0.80 
1.05 
1.40 
1.80 
2.50 
4.00 



9.14 
1.20 
9.05 
1.12 
1.42 
1.12 
7.66 
6.39 



xlQ-^ 
xlO-i 
xlQ-^ 
xlO-i 
xlQ-i 
xlO-i 
xlO"^ 
xlO"^ 



1.57 
1.13 
1.31 
1.49 
1.73 
2.88 
3.53 
8.55 



xlQ-^ 
xlO"^ 
xlO"^ 
xlO-^ 
xlQ-^ 
xlO-^ 
xlO"^ 
xlO"^ 



2.10 
1.49 
1.63 
1.25 
1.07 
1.19 
1.20 
1.03 
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